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Abstract 

Boltzmann equations and their matrix valued generalisations are commonly used to 
describe nonequilibrium phenomena in cosmology. On the other hand, it is known 
that in gauge theories at high temperature processes involving many quanta, which 
naively are of higher order in the coupling, contribute to the relaxation rate at leading 
order. How does this accord with the use of single particle distribution functions in 
the kinetic equations? When can these effects be parametrised in an effective de- 
scription in terms of quasiparticles? And what is the kinematic role of their thermal 
masses? We address these questions in the framework of nonequilibrium quantum 
field theory and develop an intuitive picture in which contributions from higher order 
processes are parametrised by the widths of resonances in the plasma. In the narrow 
width limit we recover the quasiparticle picture, with the additional processes giv- 
ing rise to off-shell parts of quasiparticle propagators that appear to violate energy 
conservation. In this regime we give analytic expressions for the scalar and fermion 
nonequilibrium propagators in a medium. We compare the efficiency of decays and 
scatterings involving real quasiparticles, computed from analytic expressions for the 
relaxation rates via trilinear scalar and Yukawa interactions for all modes, to off-shell 
contributions and find that the latter can be significant even for moderate widths. 
Our results apply to various processes including thermal production of particles from 
a plasma, dissipation of fields in a medium and particle propagation in dense mat- 
ter. We discuss cosmological implications, in particular for the maximal temperature 
achieved during reheating by perturbative inflaton decay. 
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1 Introduction 



Many important features of the observable universe can be understood as the result of out 
of equilibrium processes during the early stages of its history, when it was filled with a hot 
primordial plasma. This includes the decoupling of the cosmic microwave background, the 
creation of light elements, dark matter production, baryogenesis, and, in inflationary cos- 
mology, the production of particles altogether during reheating. In many cases the relevant 
energy scales by far exceed those that can be realised in any human made experiments and 
provide an excellent laboratory to test particle physics models beyond the Standard Model 
(SM). Nonequilibrium dynamics also play a crucial role in the understanding of signals 
from heavy ion colliders. Thus, a quantitative understanding of nonequilibrium processes 
is crucial for cosmology as well as particle physics. 

Boltzmann equations and their matrix valued generalisations can accurately describe 
many nonequilibrium phenomena. They are known to suffer from uncertainties when the 
coherence lengths are large, but are usually believed to be accurate in absence of such 
effects. On the other hand, it is known that in gauge theories at high temperature re- 
summations are necessary because processes involving many quanta, which naively are of 
higher order in the coupling, contribute to the relaxation rate at leading order. How does 
this accord with the use of single particle distribution functions in the kinetic equations? 
Medium effects are often included in effective kinetic equations by the use of thermal 
masses for quasiparticles. What are the limits of the validity of this procedure? To which 
degree can the quasiparticles be treated as particles? Do the thermal masses act as kinetic 
masses? 

All these problems are related to the assumption that the system should be described 
in terms of the properties of individual particles as asymptotic states. The definition of 
these is, however, ambiguous in a dense plasma. We show that the above questions can be 
answered in a consistent and intuitive way when swapping the single particle phase space 
distribution functions as dynamical variables for correlation functions of quantum fields. 
This formalism in addition allows for a full quantum treatment of coherent oscillations and 
quantum memory effects. 

In this work, we study the relaxation of scalar and fermionic quantum fields in a large 
thermal bath. Depending on the initial state of the out-of-equilibrium fields, they either 
gain energy from or dissipate it into the bath. This resembles a large number of interesting 
phenomena including thermal production of particles from a plasma, propagation in dense 
matter and cosmological freezeout processes. It is, apart from reporting a number of new 
results summarised below, one of the objectives of this article to explain and promote the 
formalism employed here to a wider audience. It provides powerful tools to treat nonequi- 
librium quantum systems in terms of quantities that have a clear physical interpretation, 
without semiclassical (on-shell) approximations or a gradient expansion. We aim to make 
the connection between nonequilibrium quantum field theory and effective kinetic equa- 
tions transparent also to readers without background knowledge in the former. Therefore 
we try to avoid technicalities and restrain to language commonly used in particle physics. 
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Outline of this Article 



In section [2] we introduce our notation and review tlie formalism we use, including the 
physical interpretation of the various quantities. 

Section [3] is devoted to the relaxation of scalar fields in a thermal bath. In section 
13.11 we study in detail the kinematics of the relaxation via a trilinear interaction with a 
bath of other scalars. In the first part we recall the interpretation of the known result 
in the quasiparticle approximation. In the second part we derive a formula that includes 
corrections beyond this approximation, using resummed perturbation theory. These appear 
as off-shell contributions and give rise to apparent violation of energy conservation in the 
quasiparticle picture. In section 13.21 we consider the case that the scalar is coupled to a 
bath of fermions with gauge interactions via a Yukawa coupling. In section 13.31 we discuss 
the kinematic differences between relaxation via 3- vertices and 4- vertices, using the simple 
example of a quartic interaction. 

In section m we compare the contributions to the relaxation rate from decays and scat- 
terings of real quasiparticles to off-shell contributions. 

Section [5] is devoted to the relaxation of a fermion with Yukawa coupling. In 15.11 we 
give exact expressions for the nonequilibrium two-point functions of a fermion in a thermal 
bath. In 15.21 we compute an analytic expression for the rate of relaxation via a Yukawa 
coupling in the quasiparticle regime. Details of the calculation are given in appendix IA.2I 

In section [6] we compare the time evolution of the energy density in the Boltzmann and 
field theoretical approaches. 

In section [7] we apply our results to reheating of the universe via perturbative infiaton 
decay. We discuss the possibility of an upper bound on the temperature due to closure of 
the phase space for decays by large thermal masses. 

Further implications and possible extensions are discussed in section [HI In appendix 
lA.ll we review some basic ingredients of nonequilibrium field theory used in the analysis, 
and in appendix IA.2I we derive an analytic expression for the relaxation rate of a fermion 
with Yukawa coupling in a thermal plasma that has, to the best of our knowledge, not 
been reported in the literature. 

The main new results reported in this article are the discussion of the nonequilibrium 
propagators in sections 12.21 and |6| the explicit computation of higher order contributions in 
section [3. 1.41 and the related discussion in section [3.1.3t the numerical comparison between 
leading order and resummed results in section HJ the expression for the nonequilibrium 
propagator for Dirac fermions in section |5l the analytic formula for the fermion relaxation 
rate in section 15. 2[ the comparison of the energy density in the Boltzmann and quantum 
field theoretical approaches in section |6] and the application to cosmic reheating in section 
[71 We embed these into detailed discussion that aims to draw an intuitive physical picture 
that is coherent without previous knowledge in nonequilibrium quantum field theory. 
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2 Quantum Fields and Particles 



2.1 Motivation 

Quantum Fields and kinetic Equations Many nonequilibrium processes can be 
treated in a canonical way by means of Boltzmann equations (cf. [T]) with sufficient accu- 
racy. These equations are semiclassical in the sense that they describe a system of classical 
particles which propagate freely (on-shell) between isolated interactions, the cross sections 
for which are computed from S-matrix elements in vacuum. They have been used very 
successfully to study nonequilibrium processes in a dilute, weakly coupled gas. However, 
in a dense plasma and in the presence of strong interactions the validity of these approx- 
imations is questionable. Corrections due to thermal effects have been discussed in the 
framework of Boltzmann equations [2H5]. Recently much progress has been made to derive 
consistent quantum kinetic equations from quantum field theory [5]-|24j. Most of this work 
aims to include quantum interferences, coherent oscillations and non-Markovian effects. 

Here we focus on another aspect. It is well known that the properties of particles are 
modified if they propagate in a medium. This e.g. includes Debye screening of charges 
and medium related modifications in the dispersion relations, which sometimes can be 
parameterised by thermal masses. Following early work P3H2^ on the field theory side, 
cosmological implications were soon after discussed in the context of Affieck-Dine baryoge- 
nesis in [30]. Recently the topic has received increasing interest in the context of thermal 
leptogenesis [3H5] and in models where dark matter is produced resonantly due to a change 
in neutrinos dispersion relations |31) . 

An interesting observation has been made in [S^. The authors suggest that an upper 
bound on the temperature of the universe during reheating by perturbative infiaton decay 
can be inferred when the would-be decay products acquire large thermal masses. These 
increase with the plasma temperature and eventually block the phase space for further 
decays. Such bound would have far reaching consequences as the temperature in the early 
universe plays a crucial role for the abundance of thermal relics (in particular the gravitino 
problem [23]), thermal leptogenesis [M], the fate of moduli [33] and the decompaticfication 
of extra dimensions [2S]- Indeed the mechanism would only be relevant during reheating, 
but whenever a significant amount of entropy is released by the dissipation of a component 
that redshifts like non-relativistic matter (e.g. a quickly oscillating field or non-thermal 
relic) at temperatures much larger than the corresponding mass. 

The validity of arguments based on modifications of the phase space due to thermal 
masses was challenged in |37], where the author pointed out that medium related correc- 
tions to the widths of the resonances in the final state may have a significant effect. The 
widths of decay products are known to be relevant in other contexts. In [22] it was found 
that they are essential for the damping of flavour oscillations. A ffist principles calculation 
of the lepton asymmetry generated in the decay of heavy Majorana neutrinos in [2^12^ 
shows large deviations from the result obtained from Boltzmann equations unless the finite 
widths of leptons and Higgs fields are taken into account. However, there the role of the 
width is related to the loss quantum coherence and memory, which leads to a behaviour 
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that is local in time and allows to understand the emergence of classical behaviour in the 
quantum system. The argument brought forward in [27], in contrast, refers to kinematics. 

The Particle Concept and its Ambiguities Any conclusion based on kinematic ar- 
guments about the mass, width and energy of (quasi)particles pre-assumes that the system 
can be well described in terms of the properties of single (quasi)particles that it is composed 
of. In quantum field theory this translates into the statement that integrals over spectral 
densities are strongly dominated by pole contributions. In typical collider experiments the 
language of stable or instable particles is suitable to describe the behaviour of quantum 
fields at asymptotic times because far away from the point of an interaction, elementary 
excitations of the fields propagate like free particles. It is, however, well-known that there 
exist various physical systems for which a description as a collection of (real) particles is 
not suitable. The most obvious example is a coherent classical (e.g. electromagnetic) field. 
It has furthermore been known for long |39] that the particle concept is ambiguous in a 
time dependent background, which may be provided by an external field or by a gravita- 
tional background jlQ]. In nonequilibrium systems it is also not clear that a description in 
terms of (quasi)particles is suitable in the presence of strong interactions because the def- 
inition of asymptotic states is ambiguous due to the omnipresent background plasma. For 
instance, the experimental results from heavy ion collisions have revealed that the QCD 
plasma near the temperature of hadronisation cannot be well-described by an effective 
quasiparticle model [55] . 

A consistent description of nonequilibrium quantum fields without reference to particle 
numbers or asymptotic states is always possible in terms of quantum mechanical correlation 
functions. Therefore correlation functions and expectation values of operators form the 
suitable language to treat nonequilibrium quantum fields in the early universe. Our analysis 
is based on these techniques, mainly using the formalism reviewed in [H] and applied to a 
similar problem in |6]. We will use the notation and several results of that work, which are 
summarised in section [2^ and appendix lA.il However, we believe that most of the following 
arguments can be understood without prior knowledge of nonequilibrium quantum field 
theory. 

2.2 Elements of Nonequilibrium Quantum Field Theory 

In this section we recall the ingredients of nonequilibrium quantum field theory used in 
this work, using the example of a real scalar field 0, and introduce our notation. Some 
additional formulae are summarised in appendix lA.ll For a detailed introduction we refer 
the interested reader to the references [?T] - H5] . 

We are interested in systems in which few nonequilibrium degrees of freedom are weakly 
coupled to a large thermal bath. This situation is realised in many interesting physical 
systems, including freezeout or thermal production of particles in the early universe and the 
propagation of energetic particles in a plasma. Our results can be used for a full quantum 
treatment of those. 
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As the authors of many earher works, we assume that the fields in the bath thermahse 
fast on the time scale associated with the dynamics of the out of equilibrium degrees of 
freedom. Then the background plasma can be thought of as a heat bath with a well-defined 
temperature T at each moment in time. This allows to use expressions known from thermal 
(equilibrium) quantum field theory for all correlation functions of the fields that make up 
the bath. 



The relevant Correlation Functions A thermodjTiamical system is represented by a 
statistical ensemble. In quantum field theory, this does not correspond to a pure quantum 
state, but is described by a density matrix q. The expectation value for an operator A is 
given by 

{A) = i:i{qA), (2.1) 

where we have adopted the usual normalisation Iiq = 1. The averaging (...) defined in 
(12. ip includes statistical (ensemble) as well as quantum averages. Direct computation of 
the time evolution of g is difficult in practice and only possible in few cases. Generically 
the von Neumann (or quantum Liouville) equation which governs it can only be solved 
perturbatively for a reduced density matrix with an effective Hamiltonian. In most practi- 
cal applications to date an number of additional assumptions (including an on-shell limit 
and gradient expansion) are made that lead to matrix valued Boltzmann equations which 
take account of coherent oscillation^. 

Instead of directly looking at g, it is equivalent to study the time evolution of all 
correlation functions of the fields. This can be seen in loose analogy to the Bogoliubov- 
Born-Green-Kirkwood-Yvon hierarchy in classical statistical mechanics. However, the de- 
scription in terms of correlation functions avoids all ambiguities related to the definition 
of asymptotic states and particle numbers since they are always well-defined. They allow 
to compute the expectation values for all observables at all times. Though in principle 
knowledge of the infinite tower of n-point functions is required to describe the infinitely 
many degrees of freedom of the density matrix, it is in practice often sufficient to study 
the time evolution of the mean field or one-point function {4>{x)) and two independent 
connected two-point functions. This is the case when the observables of interest can be 
expressed in terms of field bilinears in a controllable approximation. A common choice is 
given by the connected Wightman functions 

A>(Xi,X2) = (0(Xi)0(x2))e, (2.2) 

A<(xi,X2) = (0(x2)</'(a;i))e. (2.3) 

Here the subscript c indicates connected. Alternatively one can consider any linear com- 
bination of these, e.g. time ordered and anti time ordered propagators. A particularly 
convenient choice is given by 

A-(Xi,X2) = z([0(Xi),0(x2)])c = ^(A>(Xi,X2)- A<(X1,X2)), (2.4) 



"^See |45) for an example of a calculation which avoids most assumptions. 
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A+(X1,X2) = ^({0(Xi),0(x2)})e= ^(A>(Xi,X2)+A<(xi,X2)). (2.5) 

Here [,] and {, } denote commutator and anticommutator, respectively. In general the 
functions and A''= depend separately on the two four- vectors Xi = (ti,Xi) and X2 = 
(^2,^2). Here we consider spacially homogeneous fields which depend on ti, t2 and the 
relative spacial coordinate Xi — X2 only. Furthermore we are interested in the behaviour of 
fields that are weakly coupled to a large thermal bath with many constituent fields, to which 
we collectively refer as Xi. The interactions that keep the bath in equilibrium are much 
stronger than the coupling to (p. This is absolutely crucial because backreaction-effects are 
suppressed and one can compute self energies from Xi propagators only. Contributions from 
diagrams involving </)-propagators in loops are double-suppressed by the weaker coupling 
and the much smaller number of diagrams (as there are many more A'j-fields). It can be 
shown that in this case A~ is time translation invariant and depends on properties of the 
bath only [B]. 

The two point functions A^ have an intuitive physical interpretation. This can most 
easily be seen by looking at their Wigner transforms A~(ti — ^2) and A'^(ti,t2), i.e. the 
spacial Fourier transforms in xi — X2. Aq(ti — ^2) is proportional to the spectral density 
Pq{u) associated with the field operator cj), 

A-(ti-t2) = ^ r ^e-(*-*^W^) . (2.6) 
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Thus A~, which is also known as spectral function, characterises the spectrum of excita- 
tions in the plasma. To illuminate the physical interpretation of A"*" we define the energy 
momentum tensor for as T^^, = d^ipducj) — t]fiiyC\xi=o, where Xi are all fields other than (j). 
The contribution to the 0-energy e'^ = (T^) from the mode q can be written as 

iit) = \ {dtA.+^D K(tl,t2) + (0q(tl))(0q(t2))) (2.7) 

The (0) (0)-term in the brackets is due to the mean field while the A^-term is the contribu- 
tion from fluctuations that can be interpreted as particles. Thus, the statistical propagator 
can be viewed as a measure for the occupation number of the mode q that remains well- 
defined even in a dense plasma. It does not depend on any reference to asymptotic states 
or particle numbers. 

The functions A^ have to be found by solving the coupled integro differential equations 
given in flA.36p . flA.37p . known as Kadanoff-Baym equations ||33]. The solutions for scalars 



have been found in [6]. The spectral function is given by (12. 6p with 



Here H^ and H"^ are the retarded and advanced self energies. The expression (12. 8p is 
identical to the well known result from thermal field theory, but was found as an exact 
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solution for the nonequilibrium equation of motion f lA.36l) . The solution for the statistical 
propagator can be expressed in terms of A~, 

A+(ti,t2) = A+.,Aq(ti)Aq(t2) + A+.,A-(ti)Aq(t2) 

+ Kin (A-(tl)A-(t2) + A-(ti)A-(t2) 

+ r dt' r rff'A^ (ti - t')n+(t' - nA^it" - h). (2.9) 
Jo Jo 

Here the initial conditions have been parameterised as 

Km = A+(ti,t2)|ti=i2=o, (2.10) 
Kn = dt,Atiti,t2)\t,=t,=o = dt,A^{h,t2)\t,=t,=o, (2.11) 
A+.„ = d,A.Kitut2)\t,=t.=o (2.12) 

and the self energy 11"'' is defined in appendix lA.ll its Fourier transform can be related 
to the retarded self energy via f[02|) and f TOSj) . The equations fl2i:0|) - fl2i:2|) define 
initial conditions for an ensemble, which at initial time is described by a density matrix that 
can be written as a direct product g^j, <S) ^'bath- ^?bath equilibrium density matrix of the 

bath, which can be fully characterised by the temperature (and in more general cases all 
chemical potentials), characterises the initial state of 0. We assume that initially has 
Gaussian correlations only, in which case g^ is fully characterised by the five parameters for 
each mode, which can be chosen as f l2.10p -( l2.12p and the initial values (f)ci,m = (0q(t))|t=o 
and 0q4n = 9t((/'q(t))|t=o for the mean field. For A:^-^A'^-^-{A'^-J'^ = \, g^ corresponds to 
a pure state [H]. The memory integral in the last line of (12. 9p contains all non-Markovian 
effects. The time evolution of the mean field mode q is given by 

^Pq(u;)e-'-*, (2.13) 

where we assumed (0) = in the ground state. Equation fl2.13p is valid if the initial 
deviations 0q,in and 0q,in are not too large. It can be found from the finite temperature 
effective action |46|^ which leads to an effective Langevin equation for [7], that in the 
case of consideration here is equivalent to the Kadanoff Baym equations |6lll7]. Using the 
fact that = (n^)* for a real scalar field, ( 12. 8 p can be rewritten as 

-2Imn«(w) + 2uje 

^''^'^^ ~ (a;2 -col- Ren^(w))2 + (Imn^(w) + uef ' ^^'^ ^ 

This expression fulfils the well-known sum rule 

dupq{u) = 2tt. (2.15) 



^In |46j contributions to the effective action from four point functions where included, leading to an 
equation of motion that goes beyond the Langevin type, but equilibrium propagators where used for in 
loops, restricting the analysis to the regime where the relevant (p modes are close to equilibrium. 
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The poles of (I2.14p determine the spectrum of propagating field excitations. Their positions 
in the complex cj-plane depend on the self energy 11^, which depends on temperatur^. The 
self energy can be written as a sum of a temperature independent part, which coincides with 
the expression in vacuum, and a temperature dependent correction due to the medium. 
The former includes the same divergences known from the vacuum theory and has to be 
renormalised. The temperature dependent piece does not contain additional divergences 
as the medium does not affect the physics at short distances |12]. The renormalisation 
conditions can be imposed at T = 0, and the renormalised spectral function, expressed in 
terms of renormalised quantities, has the same shape as (12. 8 P [6]. We are not concerned 
with the details of the renormalisation here, but it is important that, as the renormalisation 
conditions are imposed at a particular temperature (typically T = 0), it is not possible to 
absorb any shifts of the poles at some other temperature by renormalisation. 

The correlation functions (12. 6 p and (12. 9 p are exact solutions of the quantum equations 
of motion flA.36p and flA.37p . without semiclassical or Markovian approximations. Here 



we are mainly interested in kinematic aspects, but they may also be used to study other 
quantum effects. If, for instance, there are several fields that carry an additional index 
(e.g. flavour), the definitions (12. 4p and (12. 5p are matrix valued and include correlations 
between fields of different flavour. Then the equations of motion (1A.36P and (IA.37P are 
matrix equations and the self energies (which in general are non-diagonal in flavour space) 
can give rise flavour oscillations. The formalism presented here allows a full quantum 
description of these. 

The Quasiparticle Approximation Let r2q be a pole of Pq{oj) with = Ref2q and 
Fq = 2Imr2q. If the inequality 

Tq < fiq (2.16) 

is fulfilled, the spectral density features a sharp peak at u = that can be interpreted 
as a quasiparticle with well-defined energy and a dispersion relation given by Q^- This 
interpretation requires, in addition to the condition (I2.16p . that the energies of all quasi- 
particles with the same quantum numbers are sufficiently well-separated {6Q > T) to be 
viewed as individual resonance^ Then one can approximate 

Tq^ ^ (2-17) 

and interpret Fq as the width of the quasiparticle. As the denominator of (r2.14p can be 
a complicated function of T and the wave vector q, there can be poles in addition to 
those that remain in the limit T — )■ 0. They can be interpreted as collective phenomena 



■^We have not made this dependence exphcit here and elsewhere because T here is not a dynamical 
variable, but an external parameter. All self energies, and consequently spectral functions, rates, dispersion 
relations etc. that appear in this work are temperature dependent unless stated differently. 

"^In [48] the spectral density has been used to define an effective number of degrees of freedom that takes 
the classical value if all resonances are separated and allows to interpolate into the regime of degenerate 
masses. 
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or plasma waves. In a hot QED plasma of photons and electrons, for instance, there are 
propagating low momentum modes of longitudinal photons known as plasmons, whose 
dispersion relation differs from that of their transverse counterparts. There is also an 
additional fermionic excitation with negative helicity over chirality ratio j2^, sometimes 
referred to as plasmino or hole. As the following discussion will show, the origin of the poles 
(dressed particle or collective) is irrelevant for the kinematic properties of the corresponding 
quasiparticle^. They are characterised by their dispersion relation f2q and width Fq, 
and these are the only relevant quantities for our purpose. Note also that the dispersion 
relation f2q is not necessarily similar to that of a free particle in vacuum or even parabolic. 
A temperature dependent effective mass can be defined via the curvature of f2q at its 
minimum as a function of |q|§ 

M = {df^^mm{%))~\ (2.18) 

The dispersion relation is always approximately parabolic in the vicinity of the minimum, 
though the minimum in some cases (e.g. fermionic holes in QED) is not at |q| = 0. For 
notational simplicity we denote effective masses by capital letters M, Mj without making 
the temperature dependence explicit, while small letters m, rrii refer to intrinsic (vacuum) 
masses. 

We will in the following refer to the situation = fiq as real or on-shell quasiparticles. It 
is understood that the "mass shell" for quasiparticles defined by Qq can have a complicated 
shape. Of course these particles are not "real" in the sense that they could leave the plasma 
as their properties are determined by the interaction with the omnipresent background. 
Consequently we refer to contributions from energies u ^ fl^ as off-shell contributions 
from virtual quasiparticles. Finally, we refer to the situation when fl2.16p is fulfilled and 
loop integrals sufficiently strongly dominated by energies ^ r2q as quasiparticle regime. 

It is important to notice that the above was derived without reference to any asymp- 
totic states. Furthermore, the nonequilibrium propagators f l2.6p and (12. 9 p are consistent 
solutions of the nonequilibrium equations of motion. They do not suffer from the known 
problems of approaches to nonequilibrium field theory that are based on the ad-hoc replace- 
ment of the Bose-Einstein or Fermi-Dirac distributions in the expressions for equilibrium 
correlation functions by some arbitrary functions. In particular, there are no "pinching 
singularities" and secular terms. 



Correlation Functions for Quasiparticles As is very weakly coupled by assump- 
tion, the spectrum of its excitations is similar to those in vacuum and can be characterised 

^Different authors use the word quasiparticle with different meanings. In some cases, it is used only for 
those resonances that correspond to screened one particle states, in contrast to collective excitations. Some 
authors also use it only if the dispersion relation is approximately parabolic. Here we do not distinguish 
these cases and refer to any resonance that is narrow in the sense of ()2.16p . thus has a well-defined dispersion 
relation, as quasiparticle. 

^Alternatively an effective mass can be defined as the momentum independent piece of RcII-'^ that 
comes from local diagrams, the plasma frequency flq=o when the quasiparticle is at rest, the minimal 
possible value of Xlq as a function of |q| or by demanding g^lcj^n,, = M^. 
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by quasiparticles corresponding to screened (f) particles. Then the integrals in fl2.6p . f l2.9p 
and (I2.13P are dominated by the regions near the poles, where one can approximate Pq(ci;) 
as 

with lim-r^o = i^q- Using Cauchy's theorem one obtain^ 



A,(y) - ^i^^e--^'^'^ (2.20) 



A+. ^ A+. 

A+(t; y) ^ (cos(2fiqt) + cos(fiqy)) e~^^' + sin(2fiqt)e-^''* 

A+- rothf^^ 
'^''^ ^cos(2fiqt) - cos(fiqy)) e-^-' + ^ ^ cos(fiqy) (e-"'^""/^ - e-r-'*) 



(2.21) 



and 



(0q(t)) ^ ^q.n^^^^e-^'^V. + 0q,„ COs(fiqt)e-^'^V^ (2.22) 



Here we have introduced centre of mass and relative time coordinates t = (ti + t2)/2, 
y = ti — t2 and used ( ]A.32p as well as lmll^{u) = Il~{u)/{2i), which follows from the 



definitions (1A.30P and (lA.Sip . It is easy to see that (I2.2ip does not coincide with the 
vacuum propagator in the limit T — > 0. This is because T is the temperature of the bath, 
not of 0, and the system may be prepared with arbitrary correlation functions for even 
if there is no thermal bath. This case was considered in |l9] in the context of decoherence. 

The solution fl2.2ip also differs from the free scalar vacuum propagator in the decoupling 
limit Fq — )■ 0, fiq — 7- CJq. At first sight it may be worrying that in this limit the correlation 
functions are not time translation invariant as the cos(2f2qt) depends on the centre of mass 
time. However, are not directly observable. The dependence of the (observable) energy 
( 12. 7p disappears in the decoupling limit, see section [61 The reason for the dependence on 
centre of mass time is that can be prepared with arbitrary initial correlations even in 
a free theory. This can be seen when computing (00) = Tr(^)00) using a basis \n) of 
eigenstates of the Hamiltonian H with eigenvalues En, 

Tr(^?0(ti)0(t2)) = 5^(r^|^?0(tl)0(^2)|^^) = 5^(n|^?|m)(m|0(tl)0(^2)|^^). 

n n,-m 

Here g is the initial density matrix that characterises the ensemble. If is a functional 
of the Hamiltonian only, g = V[B], then {n\g\m) = V{Em)Snm and only terms with n = 
m contribute. Decomposing 0(ti)0(t2) into creation and annihilation operators ap, aj. 



^These formulae are valid up to corrections of higher order in F/fi. The poles of /_b(w) do not give 
sizeable contributions. 



12 



it is easy to see that for n = m only the terms flpfl^, ctpCtp, which come with factors 
exp(±i(ti — t2)En), can contribute. Thus, the correlator depends on ti — ^2 only and is 
time translation invariant. If, on the other hand, the ensemble contains states that are not 
eigenstates of the Hamiltonian, terms with n ^ m appear and the combinations OpOp, a^pa^p 
contribute, which come with factors that depend on ^1+^2- The choice of initial correlations 
corresponds to a choice of quantum states to be included in the thermodynamical ensemble 
under consideration. If one imposes equilibrium initial correlations for (p with a temperature 
that is equal to the bath temperature, fl2.20p and fl2.2ip become time translation invariant 
I. 

Perturbative computations in nonequilibrium quantum field theory are sometimes per- 
formed by using the equilibrium propagators (lA.lSP with the Bose-Einstein distribution 
fsi^) replaced by some general distribution function f{oO(^) 7^ fsi^q), 

A- -> tp^iu) , A+ ^ (i + /K)) pq(a;) 

It is clear that such approximation can only be valid on timescales much shorter than the 
relaxation time ~ l/Fq because otherwise the correlation functions have to depend on 
the centre of mass time t. The fact that f l2.2ip is not time translation invariant in the 
limit Fq — )■ means that even in that case fl2.23p is not a consistent limill^ though it may 
in some cases lead to approximately correct results. In a more general nonequilibrium 
system (without the assumption of weak coupling to a large bath in equilibrium) the 
situation becomes even worse because also A~{ti,t2) depends on the centre of mass time 

(tl+t2)/S 



The Relaxation Rate Equations f|2.20p - fl2.22p show that correlations of the field are 
damped exponentially with a rate F^^H Knowledge of the initial conditions is lost after a 
relaxation time ~ F~^. Thus, Fq = — Imllq/fiq can be interpreted as the relaxation rate 
in agreement with the well-known result from thermal field theory ^T\. In general the 

^It has been discussed in j50| that it is impossible to prepare a system in equiUbrium (and thus construct 
a time translation invariant solution to the Kadanoff-Baym equations) by specifying the two-point functions 
only. The reason is that equilibrium is not a Gaussian state and higher n-point functions contain connected 
pieces. These generally enter the two-point function via the self energies, thus setting only the two-point 
functions to their equilibrium values without adjusting all other connected n-point functions does not 
give a time translation invariant solution. Here, however, self energies are computed from correlation 
functions of bath fields only, which are computed from an exact equilibrium density matrix. Therefore 
a time translation invariant equilibrium solution can be constructed by specification of the one and two 
point functions only as long as the bath is sufficiently large that backreaction is negligible. 

^One can of course always parametrise A+ = (1/2-1- f)p etc. if one allows the function / to 
have additional dependencies (including time and off-shell energies): one simply defines fqiuj^t) = 
A+(w;t)p-i(w;<)-l/2. 

"In [TU] a toy model was used to study numerically how the quasiparticle spectral function emerges 
dynamically in this case. 

^^The initial correlations in (j2.2ip are not damped with respect to y because we are studying an initial 
value problem and have assumed that the interaction is switched on at the initial time, i.e. there were no 
correlations between </> and the bath at earlier times. 
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relaxation rate is given by the discontinuity of the self energy, see appendix lA.ll in this 
case 



= -^(n^(^q) - n<(fici)) = -^n^(^q). (2.24) 

The identification of discontinuity and imaginary part that is required to derive fl2.17p from 
(I2.24P is specific to the case of a real scalar with real couplings, see (IA.35|) . The self energies 



are defined analogue to A^, see appendix lA.ll They allow to define the quantities 

= -2k^t(^n) r> = -^n>(fij. (2.25) 

Fq can be interpreted as the gain rate of thermal production of quasiparticles from a plasma 
and Fq as the loss rate due to the inverse processes. Their difference Fq = Fq — F^ then 
gives the total relaxation rate for the mode q through processes involving quasiparticles 
with energy u = fiq. The rates F^ fulfil the detailed balance ratio F^/F^ = exp(— 
as a result of the fact that the bath is in equilibrium, see ( lA.17p . ( IA.32I) . Note that this 



definition of the relaxation rate Fq does not involve any asymptotic states, the definition 
of which is ambiguous in the omnipresent plasma. 

The interpretation of Fq as a relaxation rate, i.e. the rate at which the mode q of 
(p exchanges energy with other modes, holds regardless of whether the initial occupation 
is below or above its equilibrium value. If the initial conditions in (12. 9p are chosen such 
that occupation numbers are below their equilibrium values, thermal production of quanta 
from the plasma drives the relaxation while in the opposite case the excess of energy 
stored in is dissipated into the batlj^. The relation fl2.24p is well-known from linear 
response theory |42] and has also been derived from first principles for systems where the 
occupation numbers are small at all relevant times [H]. Note that in our setup this 
interpretation holds irrespective of the size of the deviation of from equilibrium as long 
as the bath has sufficiently many degrees of freedom. In the following we use f l2.24p as the 
definition of Fq. 

The fact that the background plasma is in equilibrium allows to use the techniques of 
thermal (equilibrium) field theory to compute the self energies. There exist two different 
formalisms in thermal field theory, known as real and imaginary time formalism |32]- The 
latter, also known as Matsubara formalism |52| . is more commonly used. Here we employ 
the real time formalism, which is physically more transparent and can directly be integrated 
into the Schwinger-Keldysh formalism used to compute the nonequilibrium correlation 
functions for [13]. The reason is that in both formalisms correlation functions are derived 
from a generating functional for fields with time arguments on a contour in the complex 
time plane and Feynman rules for equilibrium fields are formulated analogue to those in 
the Keldysh formalism sketched in appendix lA.lt with (lA.lSP inserted into (lA.23P - (lA.26p . 



^^The fact that > e ^'^T^ = r< does not imply that (f> always loses energy because are just the 
rates. The gain and loss of energy for each mode also depends on the occupation number, c.f. (|6.ip 
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Figure 1: Feynman diagrams contributing to Fq, the lines represent (solid), Xi (dashed), 
X2 (dotted), \E'i (solid with arrow) and \i/2 (dashed with arrow), the gray circles symbolise 
self energy insertions for the Xi due to some interactions with themselves or other fields in 
the plasma. 



3 Relaxation of a Scalar 

In this section we compute the relaxation rate for scalars with various interactions. We 
focus on the effects that medium related changes in the dispersion relations of particles and 
plasma waves in the bath have on Fq. We will leave the energy u free wherever possible 
for the sake of generality, but we are mainly interested in the case w = f2q ~ Wq, where 
the second equality expresses that corrections to the dispersion relations of 0-resonances 
are small due to the weak coupling (while they have to be taken account of for all other 
fields in the bath that is kept in equilibrium by stronger interactions). To be specific we 
chose the Lagrangian 

C = ^9^09^0 - im202 ^ ^ f^^,x^^^x^ - \m\x\ + (^^ " m.) 

-#XlX2 - Y<^ (^1^2 + ^2^l) + {-\'^^''^ + '^Xdnt + /:^,int) • (3.1) 

Here ^ is the out-of-equilibrium field. The Xi scalar and the fermionic fields that 
are part of the thermal bath. The bath is thought to contain many more fields Xi to 
which the Xi ^-^id couple via interactions contained in the terms C^r-^^x,^ -^^iint- Their 
particular form is of no relevance at this point, they are only assumed to be strong enough 
to restore equilibrium in the bath fast on the time scale associated to the dynamics of 
0. Backreaction can be neglected if the bath is sufficiently large. The Lagrangian (13. ip 
may be regarded as a toy model, but it allows to study all relevant kinematic effects while 
avoiding the complications of gauge theories at finite temperature. 



3.1 Trilinear Scalar Interaction 

We first consider the trilinear coupling g(f)XiX2 in (13. ip . The self energy is then given by 
the diagram shown in figure [1^). We use the fact the n~ can be expressed in terms of 
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n^, which in the real time formahsm of thermal field theory at leading order is given by a 
single diagram (combination of contour indices). Imll^ to leading order is given by 

ImnJ(a;) = = ^/^^(a;)n<(a;) 

4 

■0 - w 



-y/B'(^) / ^Afp(Po)A<p_,(a;-|>o 



= - Y y (1 + /(Po) + /(t^ - Po)) Pip(po)p2p-q(w - Po)- (3.2) 

Here is the Bose-Einstein distribution /b(i^) = (e^"^ — 1)~^ and are the spacial 
Fourier transforms of the connected thermal Wightman functions {Xi{^2)Xi{^i)) c defined 
analogue to (12. 2p . By means of flA.18|) it can be written as 



Afp(po) = /b(Po)p»p(Po)- (3.3) 

We have also used the identity fB{ci)fBip — (^){fB{b))~^ = 1 + /^(a) + /^(fe — a). The 
spectral densities Pip(po) for the fields Xi ^-re defined analogue to (12. 8p . In the remainder 
of this section quantities that carry an index j (such as rrii, Mi, pip{po) ■ ■ ■) are always 
associated with the bath fields Xi while quantities without such index (m, M, Pq{u}) ■ ■ ■) 
refer to (p. Equation (13. 2p shows that the quantities that determine Fq are the spectral 
densities of the bath fields Xi- 

3.1.1 Relaxation Rate at leading Order 

A leading order result for Fq can be obtained by inserting free spectral densities, 

pf-(po) = 27rsign(po)5(p' - m^), (3.4) 

into fl3.2p . Equation (13. 4p can easily be found from fl2.8p in the decoupling limit. Then 
reads 

Imn^(a;) = ^ ' 



' 2 J (27r)Mw2Wi 

X (((/l + 1) (/2 + 1) - /1/2) (<^(^ -^1- ^2) - 5{UJ + UJI + ^2)) 

+ ((/i + l)/2 - (/2 + l)/i) -00, + U2) - S{u + u;^-U2)y 



(3.5) 



Here ooi = (p^ + miy^'^, 002 = ((q — p)^ + miy^'^ and /j = fsi'^i)- The well-known 
expression (13. 5p has a clear physical interpretation [51]. The first ^-function in the first 
line corresponds to decays and inverse decays -H- X1X2, see figure [2|l). Looking at the 
statistical factor, one can easily confirm that the detailed balance relation is fulfilled, a 
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Figure 2: Diagrams obtained from cuts through the diagram in figured^). 

consequence of the fact that the bath is in equihbrium. This factor leads to an amphfi- 
cation of the rate compared to the vacuum value due to induced transitions, typical for 
bosons. The 5-function implies energy conservatioij^. The second (5-function in the first 
line corresponds to the creation of real Xi ^^nd a virtual quantum with negative energy 
from the vacuum, see figure [2)3). It comes with the same statistical factor and conserves 
energy as well, but never contributes in the interesting case u = Qq- The second line 
corresponds to processes of the type Xi ^ emission or absorbtion of a </> quantum 

by the bath, see figure [2t) and d). Again, detailed balance as well as energy conservation 
remain valid. This contribution involves interactions with real quanta from the bath and 
leads to a relaxation mechanism similar to Landau dampin^^. In the vacuum limit T — )■ 
the second line vanishes, and the first line reduces to the known decay rate in vacuum. 
The integral in (13. 5p can be solved analytically [3153] and gives 

ImU^iu) = ao{q) + aa{q,T) + a,{q,T). (3.6) 

Here ctq is the contribution due to the decay process (p — )■ X1X2, 



^o{q) = - -77, — ^sign(u;)6'(g^ - (mi + mg^^^ 



X 



{{q'f - 1q\m\ + m\) + (m? - m^)^) ^ , (3.7) 



aa{q,T) is an additional temperature dependent contribution from such processes that 
includes induced transitions and inverse decays. 



2 



^a{q, T)=- J sign(a;)g(g^ - (mi + ms)^) 

^ f^^ + ("^1 ^ ^2)) , (3. 



^•^The quasiparticle momentum is conserved exactly at all orders because we restricted the analysis to 
systems that are invariant under spacial translations. 

^''Traditionally the term Landau damping refers to the interaction of a particle with a wave or mean field. 
In the following we use this term more generally for any process in which </> interacts with the background, 
i.e. that has quanta of fields in the bath (other than (p) in the initial state. 
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and (Tfe(g) is the contribution from Landau damping 



(g, T) = - ^ sign(u;)g((mi - ms)^ - q 



1 _ e~'3|^^ 



+ (mi O ms) . (3.9) 



We have used the abbreviations 



oo± = |^(g' + - ml) ± ((g^ + _ ^2^2 _ 4^2^^, | ^ ^g^^Q^ 

where q is the four vector q = (a;,q). Kinematic restrictions known from vacuum remain 
vahd, energy and momentum are conserved in individual processes. The differences to 
the vacuum he in the statistical factors and the fact that additional processes occur. The 
plasma can provide real Xi quanta in the initial state that make inverse decays and Landau 
damping possible. 



3.1.2 Resummation 

The energy conserving 5-functions originate from the use of free spectral densities (13. 4p . 
They define submanifolds in the integration volume of (13.21) in which the factors Afp(po) 
and A^p_q(a; — po) are non-zero, the mass shells of Xi particles. Only if those intersect at 
some point in the integration volume, the integrand is non-zero. In vacuum this is well- 
known from the optical theorem and the cutting rules it implies |54| . These imply that 
n~ is non-zero only if all particles in the loop integral can be on-shell somewhere in the 
integration volume. The finite temperature generalisation of these rules |55|^ loosely state 
the following: Cut the diagram in all possible ways and consider all possible amplitudes 
that can be built from the pieces by interpreting their external legs as initial and final 
state particles. Cut propagators are replaced by or A^, depending on the momentum 
flow (for fermions S^). The spectral densities appearing in A^, see (|A.18l) . ensure that the 



amplitudes only contribute if all external particles can be on- shelff while the statistical 
factors lead to factors J'b for initial and 1 + /s for final state particles. The former ensures 
that in the limit T — contributions from processes with particles provided by the heat 
bath in the initial state vanish, the latter reflects the effect of induced transitions for bosons 
(for fermions there is a suppressing factor 1 — fp, see (lA.lOp ). Therefore, as in vacuum, 
only those cuts through the original diagram contribute for which all cut propagators are 
on-shell simultaneously somewhere in the integration volume. The additional contributions 
have an intuitive physical interpretation as processes in which the external particle of the 
original diagram engages in reactions with real quanta from the plasma. 



^'"'The original formulation of the finite temperature circling rules |57| did not make this interpretation 
in terms of cuts and products of amplitudes obvious. 

^^So called pinching singularities in the single amplitudes, see e.g. |56) . are regularised as we in the 
following argue that a consistent treatment requires resummation of all propagators. 
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Figure 3: An example for a ladder diagram, the lines represent (solid), Xi (dashed), X2 
(dotted, small spacing) and some other unspecified field X (wavy line) in the bath. 



Equation fl3.5p only accounts for the effect of the plasma on </>, but not on the Xi- 
This is inconsistent as we assumed that interactions C^^i^t are much stronger than the 
coupling of the Xi to 0. Their effect can be included by using resummed propagators in 
(13. 2p . which corresponds to inserting dressed spectral functions of the form (I2.8|i ^1. Note 
that this procedure, though being a consistent resummation of the propagators, does not 
take into account all possible contributions to the self energy. It e.g. neglects final state 
interactions coming from ladder diagrams as shown in figure El In the following focus 
on the case that diagrams of the type shown in figure [Tb) dominate the self energy in 
some controllable approximation. The quantitative validity of this approach depends on 
the details of the interactions that keep the bath in equilibrium. It is e.g. justified when 
there are either many more or stronger interactions giving rise to diagrams of this type. 
We believe that most of our conclusions qualitatively hold when this restriction is liftec^j, 
but computations become more cumbersome and a widely analytic treatment as we will 
perform in the following in general is not possiblj^ 



Quasiparticle Regime The dressed spectral densities Pip{po) for the Xi can be compli- 
cated in general as we have not specified their interactions. Let us for a moment assume 
that all poles fulfil (I2.16p . Then the Pip(po) are characterised by a number of narrow peaks 
at energies uj = if^fp, corresponding to quasiparticles with dispersion relation fifp and 
width rfp. Here the index i indicates that we refer to a pole of the spectral density Pip{uj) 
of Xi and the index I labels the various poles. 

Near each pole the spectral density can be approximated by a Breit Wigner function 
in the usual way. 

In the limit of vanishing widths Ffp, this converges to a sum of 5-functions of the form 

^^The need to employ resummed propagators to obtain consistent results is common in finite temperature 
field theory [4l[4ll|421[58H60] . 

^^This is e.g. confirmed by the results found in |4] while this work was in progress. 

""^^If one employs nPI techniques as e.g. discussed in |41) all contributions at a given order are automat- 
ically taken into account. Here we perform the resummation "by hand" in order to make maximal use of 
the simplifications due to the weak coupling to a large thermal bath. A more general, but less transparent 
resummation scheme has also been proposed in ??. 
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~ (5(a;^ — (f2fp)^). Then the previous argument can be repeated, the 5-functions define 
submanifolds in the integration volume of (13. 2p on which the integrand is non-zero. Only if 
these submanifolds intersect somewhere in the integration volume, there is a contribution to 
the integral. This can be interpreted as energy and momentum conservation in scatterings 
and decays of real quasiparticles. Note that, as pointed out previously, this argument is 
independent on whether the quasiparticles are screened particles or collective excitations. 

In the simplest case, when the only poles of Pip{po) are ±(p^ + Mf Y^"^, Fq is given by 
(I3.6p - (l3.9p with intrinsic masses replaced by thermal masses, rrii — Mj. This allows to 
define two critical temperatures and as the solutions to 

= (M - Ml - Ma)^.^^^ (3.11) 
= {M - \Mi - M2\)T=f, . (3.12) 

For T < Tc, (13. 7p and ( 13. 8 P contribute to Fq (decays and inverse decays O XiX2)- We 
refer to this as case (a). For T > T^, Fq is given by (13. 9p due to processes (pXi ^ Xjj 
case (b). For < T < none of these processes is allowed and Fq vanishes in this 
approximation, case (c). The generalisation of the definitions Tc and Tc and the cases (a), 
(b) and (c) to situations with more complicated dispersion relations is straightforward and 
we will in the following use them in this generalised sense. 

When inserting the full spectral densities including continuum parts and finite widths, 
the integrand in (13. 2p is non-zero everywhere in the integration volume. However, as long 
as ( I2.16P remains fulfilled for all poles, F^p <C ^Li, ^he integral is strongly dominated by 
the regions in which all quasiparticles are on- shelH The region in the integration volume 
where Pip(po) peaks extends only a small distance ~ Ffp away from the quasiparticle mass 
shells Po = (^ip)^ Contributions from regions where one or more quasiparticles are off- 
shell are suppressed, the suppression is stronger in regions where more quasiparticles are off- 
shell. Hence, the integral (13. 2p only receives a large contribution if these quasiparticle mass- 
shells intersect or get very close to each other. Otherwise it is non-zero, but only receives 
contributions from regions in the phase space where at least one of the quasiparticles in the 
loop is off-shel|^. The suppression is analogue to the suppression of processes in vacuum 
that involve virtual intermediate states which are much heavier than the energy of the 
process. 

In the case discussed here, the relaxation of a single real scalar field, there is only 
one relaxation time scale ~ 1/Fq per mode. In more complicated systems there can be 
several time scales related to kinetic equilibration of different species, effective chemical 
potentials or quantum coherences (i.e. correlations between different fields). In [T7] it was 
found that the on-shell approximation can lead to "spuriously conserved quantities", i.e. 
quantities that appear to be conserved though there is no corresponding symmetry in the 
Lagrangian. This is analogue to the apparent lack of relaxation in regime (c) observed in 
the on-shell approximation here. In jl^ it was pointed out that the "spuriously conserved 

^°Note that the loop integral in this case is not dominated by momenta |p| ~ T, but by the on-shell 
regions. 

^^If there are several collective resonances, the line in the loop can represent a propagator for any of 
these. The decay into holes has e.g. been studied in [5]. 
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quantities" disappear when the on-shell approximation is hfted. They do, however, relax on 
time scales much longer than the kinetic equilibration because the processes that violate the 
conservation involve off-shell quasiparticles. This is in complete analogy to the suppressed 
relaxation rate in regime (c) in our case. 

The fact that the dispersion relations change with T implies the phase space for various 
processes is dynamical and can open or close at critical temperatures. (13. 6p shows that 
this may happen rather abruptly leading to the "thresholds" along the temperature axis 
that are visible in figure [HI This can lead to interesting dynamics as the temperature itself 
can be affected by the relaxation processes (see e.g. the example discussed in section [7]). 

Broad Resonances If, in contrast, the widths of the Xi are large, there can be a sig- 
nificant overlap of the spectral densities in (13. 2p even far away from the on-shell regions. 
Figure [8] shows that in this case the thresholds get smeared out. Fq changes rather con- 
tinuously with temperature due to the smooth shape of the broad resonances. There is no 
critical temperature at which Fq changes abruptly. 

A resonance that has a narrow width at temperatures well below its mass may become 
broad and loose its identity at hight T. This phenomenon, known as melting of a peak, is 
theoretically well-studied and has been observed experimentally [65] for mesons. Quantita- 
tive computations in the high temperature regime are difficult due to the poor convergence 
of the perturbative series even in weakly coupled theories. They usually rely on resum- 
mation methods, lattice computations or effective field theories of lower dimensionality. 
Recently, also the use of gravity-field theory dualities has been explored. 

Physically the broadening of meson resonances is of course related to the approach to 
the QCD crossover and the dissociation of the meson states. In theories of fundamental 
particles, where broadening can also occur, it means that the quasiparticle cannot be 
viewed as an entity that is well-separated from its environment. It is taken off-shell by the 
statistical fluctuations of its energy due to the interactions with the dense medium. The 
apparent energy non-conservation in the quasiparticle decay is not surprising as only the 
combined energy of system and environment is conserved. From a quantum viewpoint this 
means that processes generally involve many quanta. From the previous considerations 
and figure |H] it is obvious that also the opposite can happen: The width of (p suddenly 
becomes more narrow at T^. 

3.1.3 Analytic Structure of the Spectral Density 

The different kinematical constraints can be studied in terms of the analytic structure of 
the spectral density as a function of squared four momentum, schematically shown in figure 

s 

Leading Order In vacuum Pqioj) has a continuous part above the lowest multiparticle 
threshold g^^^ = (mi -|- 7712)^. In this region it inherits the discontinuity of the self en- 
ergy across the real axis given by (13. 7p . This discontinuity can be computed by cutting 
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(b) (c) (a) 

Figure 4: Analytic structure of the spectral density of as a function of for a self energy 
from trilinear interactions (i) in vacuum, (ii) at finite temperature with bare propagators 
and (iii) at finite temperature after resummation. The zigzag lines represent discontinuities 
and the dots the positions of poles in the cases (a) > q'fi^i, (b) < q^^^ ^'^^ (c) 
^th2 < < Qthiy where Q = (fiq, q) is the on-shell energy-momentum four vector. 



the diagram in figure [T^). Poles that lie below qthi lead to singular, 5-function shaped 
contributions that are interpreted as stable particles and possibly bound states. 

In a medium the diagram in figure [T^) also contains the processes (pXi Xj^ leading 
to an additional discontinuity below a threshold g^^^g- ^^e simplest case, when the 
effect of the medium on the Xi dispersion relation can be parametrised by momentum 
independent thermal masses, it is given by g^^^g = (^i~^2)^; see ( I3.9|l l^. The discontinuity 
dresses any poles below qth2 with a finite width. Defining the four vector Q = (fiq, q) one 
can distinguish the cases (a), (b) and (c). In case (a) the position of the pole f2q is 
such that > q^^^ and relaxation is driven by decay, see (13.71) and (13. 8p . In case (b) 

< qfi^2 dissipation is driven by Landau damping, see (13. 9p . In between, there is a region 
^th2 < < Qthi which none of these processes is possible for real quasiparticles, case (c). 
A similar structure is found in solid states, where it determines the conductivity. A pole in 
region (c) leads to quasiparticles that can move freely over large distances. It depends on 
T which of these cases is realised. The critical temperatures Tc and mark the transitions 
from (a) to (c) and from (c) to (b), respectively. The three temperature regimes are clearly 
visible in figure [H 

Multiparticle Final States and Scatterings Cutting the resummed diagram in figure 
[Ij)) through the Xi self energy insertions leaves pieces that correspond to the decay of (p 
into multiparticle final states that are composed of whatever particles Xi appear in those 
insertions, and the inverse processes, see figure [5^). The Xi appear as intermediate state. 
These decays contribute if they are kinematically allowed, but they are suppressed by 

^^Note the difference to [5T], where this discontinuity was found to vanish for < — \m\ — m||. 
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Figure 5: An example for diagrams obtained by cutting the diagram in figure [TJo) if X2 
couples to an unspecified field X in the bath, lines represent (solid), X\ (dashed), Xi 
(dotted) and X (wavy) propagators. 



additional factors of the involved coupling constant. In the case (c), when the intermediate 
Xi cannot be on-shell, they are also suppressed by the smallness of the Xi width compared 
to its energy, which corresponds exactly to the off-shell suppression previously discussed. 

In addition, the diagram in figure [TJo) also contains scattering processes with Xi in the 
initial and final state, see figure ISjo). In contrast to the decays, they are kinematically 
always possible for some momenta, thus there is always a discontinuity along the entire 
real axis. 

We should add that the diagram in figure[Ifc) contains scatterings of the type Xi0 ^ Xj'P- 
These processes do not change the "number" of </> quasiparticles, but exchange energy 
between </> and the bath. We have previously omitted this diagram (and will do so in the 
following) as it is a "backreaction" term that is suppressed by two additional powers of the 
weak coupling g and the large number of degrees of freedom in the bath. 

In general there are also ladder diagrams as e.g. shown in figure [3l which involve other 
fields than 0. They contain vertex corrections as well as contributions from scatterings 
and are more important than the diagram shown in figure [It) because the interactions in 
the bath are stronger than g and there may be various couplings that give rise to diagrams 
of this type. Here we consider the case when contributions of the type shown in figure [Do) 
dominate the self energy, as discussed at the end of the first part of section 13.1.21 This 
allows to parametrise the effects of all higher order processes involving more quanta in the 
initial and final state in the quasiparticle widths of the x@- The fact that these processes 
can be included without reference to asymptotic states in the omnipresent plasma is one 
of the benefits of the formalism employed here. 

Multiple Scatterings All diagrams that contribute to scatterings contain additional 
vertices. Furthermore, in case (c) at least one of the intermediate propagators must be off- 
shell, leading an additional suppression of Fq between Tc < T < T^. At low temperatures 

^•^The off-shell suppression observed in |18j that leads to a delayed relaxation of quantities which spuri- 
ously conserved in the on-shell limit can also be interpreted in this manner: These quantities relax only 
via scatterings or decays with multiple quanta in the final state, which can be extracted from cuts through 
the dressed Xi propagators. 
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scatterings are subdominant with respect to the decays. At at high temperature the small 
coupling suppression is compensated by the large occupation number^ and the scatterings 
become increasingly important. Yet, there is another subtlety. The suppression of diagrams 
involving more vertices can be compensated even in region (c) if the involved quanta 
are collinear and almost on-shell. Interferences between subsequent scatterings are not 
negligible, similar to the Landau- Pomeranchuk-Migdal effect [6^, see also |6T1463] . This 
becomes increasingly important when the temperature is larger than all involved masses 
and therefore affects the regions (b) and (c). Physically it means that the relaxation is not 
driven by processes involving only a handful of quasiparticles, but a large number of quanta 
from the bath. Intuitively it is clear that multiple scatterings should become more frequent 
at high temperature: even a small coupling leads to a large scattering probability if the 
density of scattering partners becomes sufficiently high. In this regime our approximation 
that the self energy is dominated by diagrams of the type shown in figure [Do) generally 
becomes increasingly bad above some temperature even if it is well-justified at low T. 
A consistent scheme to resum contributions from processes involving many soft collinear 
quanta has been presented in [6l]- While this work was in progress it was applied in [1] 
and found that these effects can not only lead to an increase of the rate in region (b) by 
a factor of order one, but also almost completely overcome the suppression in region (c). 
This is in agreement with our results shown in figure [HI 

The increasing importance of scatterings with temperature is of course known in the 
Boltzmann approach, but we would like to emphasise that the field theoretical consider- 
ations show that for a consistent treatment at temperatures larger than the masses, the 
use of resummed perturbation theory is necessary. This corresponds to a summation of 
infinitely many diagrams with large numbers of particles in the initial and final state in 
the Boltzmann approach |15| . 

It follows from the above that there can be no stable excitations in a plasma. Of 
course this is expected because it is clear that a particle propagating in a medium can 
always engage in scattering processes. If this was not the case, one could easily chose the 
parameters of the Lagrangian f l3.ip in a way that (p in some temperature regime does not 
exchange energy with its environment and relax despite their coupling and the possibility 
of scatterings. However, in case (c) the exchanged Xi is off-shell and the rate for such 
processes suppressed in the quasiparticle regime, thus can be long lived and stable on 
the relevant timescale. In our setup, where ladder diagrams are neglected, all higher order 
contributions are encoded in the energy and momentum dependent finite widths of the 
Xi- Therefore the deviation from the leading order contribution from decays and inverse 
decays of single quasiparticles can be parameterised in terms of Fj/f^j. 

3.1.4 Off-Shell Contributions 

In the following we concretise the above considerations and derive an explicit formula for the 
relaxation rate Fq for the zero mode of (p including off-shell contributions. For |q| = 0, the 

^^This is a physical reason for the well-known problem of poor convergence of the perturbative expansion 
in the coupling at high temperatures. 
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only spacial momentum in (13. 2p is p and we can drop the momentum index. For simplicity 
we restrict the analysis to the case that the only poles of the spectral densities for the Xi 
are those that correspond to the screened one particle states. We furthermore assume that 
the dispersion relations are in good approximation parabolic and define Cli = Qi + ^Ti with 
Qi = (p2 + M2)V2^ where Mi are the thermal masses for the fields Xi- The simplest way to 
realise this is to assume that the leading order contribution to the self energies Ilf of the 
fields Xi comes from local diagrams as e.g. shown in figure [7^) that lead to a momentum 
and energy independent mass shift. We take Fj 7^ 0, but Fj ^ Qi. With 03.21) and fl2.8p 
one can then write 



Imn^(c.) = -y / (1 + /(Po) + - Po)) 

2Imnf(po) 

2Imnf (w - Po) 



(3.13) 



We have dropped the e prescription as we know from the previous considerations that at 
finite temperature Imllf^ is non-zero along the entire real energy axis. The Pq integral can 
be evaluated using Cauchy's theorem. The integral is dominated by the on-shell regions. 
As the smooth functions Ilf do not change much within the narrow peak regions along the 
Po axis and Imllf in the denominator is irrelevant everywhere else, one can replace po by 
its pole value in the argument of Imllf . However, one cannot directly replace the spectral 
densities using (I2.19P before integration as this replacement is only valid if the function that 
Pi is multiplied with under the integral has no additional poles. Here one has to proceed 
pole by pole and ImH^ is evaluated with different arguments. It is straightforward to 
obtain 



f d'p f[l + fB{^i) + fB{uJ-ni))lmnf{ni)lmU§{uj-n 



2 J (27r)3 \^ ^^^^^^ l^j^^^ _ ^^^2 -n^y + (lmU^{u - l^i))'^ 
fBiu + ni) - fBi^l)) ImHf (-fit)Imn^(u; + n*,) 



fBiuj + (12) - fBi^2)) IniHf (w + fi2)Imnf (-^2) 



fi2r2f^2 [{CO + ^2)2 -niy + (imu^iu + ^2)) 



1 + fB{ni) + fB{uj - n*) imHf (w - n*)imn§{n*) 

^- 7 — 7^ 1 • (3.14) 

Q2T2^2 1 ((w - n*2y - nt) + (imnf (^ - n*S 
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Here we have used the identity z) = —1 — /^(-z). The angular integrations can be 
performed straight away. We now expand the numerators in the small widths Fj. With the 
identity /^(x) = — + and replacements of the form Imllf — — fiiFi 

and Imllf — ?■ — Fi one can simplify the expression to 



Imn^H = r f d\p\ / Imn^(a; - il,) (1 + /b(1^i) + fBjco - il,)) 
2 y (27r)2 |^j^(^_^^)2_f^2j% (^Ij^n^(^-fii))'^ 

Imn^(a; + fii) (/^(a; + fii) - /^(fii)) 



^2 (^((w + ^2)2 - 1^?)' + (imnf (a; + (12) 
Imnf (a; - ^2) (1 + /b(^^2) + - ^^2)) 



fi* 1^(^(0; - fi^)2 - fif) V (^Imnf (w - fi ' ^ 



(3.15) 



Each line in fl3.15p is proportional to a factor ImFlf /fij ~ Fjfij/fijCj and is thus suppressed 
by the smallness of this ratio in comparison to the characteristic quasiparticle energy fij. 
This suppression can be cancelled if the denominator is very small for some values of p. 
One can easily see that this happens if either u ^ ±(fii + fi2) (first and fourth line, 
corresponds to case (a) ) or w ~ ±(fii — fi2) (second and third line, corresponds to case 
(b) ). These conditions exactly require that processes (p o X1X2 (first and fourth line) 
and Xi ^ 4'Xj (second and third line) are possible on-shell, but with particles replaced 
by quasiparticles. Figure [8] demonstrates the efficiency of the suppression in the forbidden 
region corresponding to case (c) by a numerical evaluation of ( I3.15P for different Fj. For 
Ti/Mi ~ 10~^ it is efficient and Fq is indeed an order of magnitude smaller than the vacuum 
decay rate and two orders of magnitude smaller than at T < T^. However, already for a 
moderate width Ti/Mi ~ 10~^ the relaxation by off-shell processes can be as efficient as 
the vacuum decay because of the effect of induced transitions. 
Using the replacement 

ImHifx) . , . . 

— )■ — sign(a;j7ro(. . .) 



(...)2 + (Imn,(x))2 



it is straightforward to verify that in the limit of vanishing ImHf , f l3.15p converges to (13. 5p . 
but with the intrinsic masses m, rrti replaced by temperature dependent thermal masses 



^'"'For large |p| one has Vli « Vlj. For iioii-relativistic |p| it is possible to have VLi ^ Then one of the 
Xi quasiparticles is effectively massless in comparison to the other and there is no forbidden region (c) as 
Ithi ~ qth2 and dissipation is always possible on-shell. 
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M, Mi. Without the simphfying assumption of parabohc dispersion relations one would 
also obtain a formula containing energy conserving (5-functions. In that case these would 
depend on the full dispersion relations and could not be obtained by replacing intrinsic 
masses by thermal masses in f l3.5p . 



3.2 Yukawa Interaction 

The relaxation rate from the Yukawa coupling can be computed from the diagram in figure 
[If). Analogue to ( 13. 2p it is given by 

ImU^iu) = ^f^\uj) I ^tr {S<,iPo)S>^^^{po-u)) 

= — ^/b^(^) J ^^/F(po)(l-/Fbo-w))tr(pip(po)p2p~q(po-w)), (3.16) 

where S^{pq) = — /F(po)pjp(Po) is the spacial Fourier transform of the thermal Wightman 
function for fermions, see ( 15. ip and ( lA.lOp . Here fpiPo) = + 1)~^ is the Fermi-Dirac 
distribution and Pip(po) the fermion spectral density. 



3.2.1 Fermion Spectral Density 

Explicitly Pip(po) reads 



The pole structure of (I3.17p . which determines the spectrum of fermionic resonances, is 
complicated in general |271[29| [66]. Gauge and Yukawa interactions are known to give rise 
to additional poles that correspond to collective excitations known as holes or plasminos. 
We consider a fermion bath that is kept in equilibrium by an abelian gauge interaction of 
strength c@ which, for simplicity, couples in the same way to \E'i and \I^2- We furthermore 
assume that all resonances in the bath have a narrow width that can be neglected. We 
consider the two limits mi ^ aT and vcii ^ aT in which Pp(po) is well-known. 




m.j ^ aT : In this case the resummed spectral density in the quasiparticle regime with 
negligible width can be computed using the hard thermal loop approximation |59] . The 
well-known result reads [28 | B2 ] 

Pp(Po) = ^ ((70 - P7)P+ + (70 + P7)P-) (3.18) 
where p7 = Pi7i/|p| is the unit vector in p-direction and 

9±{p)^M9t\p) + 9T{p)) (3.19) 

^^The behaviour in noii-abelian gauge interactions and Yukawa couplings is similar |42II67| . the spectrum 
in the supersymmetric case has been studied in [BHj . 
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with 
and 



(3.20) 



x] 



X 



1 =F X ± ( (1 =F a;) In 



X + 1 



a: — 1 



±2 +7ry(lTa;) 



(3.21) 



Here x = Pq/\p\, = (a;T)^/(4p)^ and the residues are 

4y2p2 

The dispersion relations 0+ and — Q_ are the solutions to 



(3.22) 



= po-\p\[l + yU{l-x)\n 



x + 1 
X — 1 



(3.23) 



There are a small continuous contribution for Pq < |p| from p'^^^{p) and four poles in 
p^°''^(p). This complicated structure makes it obvious that simply replacing bare masses 
by thermal masses for fermions is not a valid procedure even in the quasiparticle regime. 
However, the poles can be interpreted as quasiparticles with dispersion relations determined 
by the p-dependence of their position. 0+ is interpreted as dressed particle and D._ as 
hole or plasminc0. The pole contributions dominate the integrand in (I3.16p . which again 
is given by a product of spectral densities. Thus, the entire discussion from section [3?T] can 
be repeated. The energies of quasiparticles are given by 0±, and energy is approximately 
conserved in reactions between them. 

If quanta are energetic enough, they can decay or be created by inverse decays, 
(j) -f-)- If there is a combination of two fermionic quasiparticles and (f) such that one 

can decay into and the other, \E^j ^ ^^'jEfl, Landau damping is at work. In absence of 
any combination of and \i/2 excitations such that one of them decay into the other 
on-shell, ( I3.16P only receives contributions involving the continuous parts of the spectral 
functions and Fq is suppressed as in case (c) in section 13.11 If large widths are involved, 
the kinematic thresholds get smeared out as shown in figure [8] for scalars. 

The situation simplifies for large momenta |p| ^ aT, meaning <C 1. Then one can 
approximate 



y2(l + 21n(y)) 



a ^y ^exp(-i/ ^) 



(3.24) 



''Note that the hole has a negative chirahty over helicity ratio even for fermions that are massless in 



vacuum. 
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\1/ can symbohse any fermionic quasiparticle, screened particles or holes. The decay into holes has e.g. 
been studied in [S]. 
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and 

a+ ~ |p|(l + 2y2) ~ (p2 + 97^2)1/2 a_ ~ |p|(l + 2a-iexp(-?/-2)) (3.25) 

with DJlf = aT/2. Thus, for |p| ^ aT the dispersion relation 0+ for the screened particle 
approaches that of a particle with mass ~ 931/ = aT/2. This thermal mass only appears 
in the dispersion relation and does not break the chiral symmetry [28], as one can see from 
f l3.18p . DJtf is sometimes referred to as asymptotic mass and differs by a factor a/2 from 
the plasma frequency at |p| = 0. Equation (I3.25P shows that the holes effectively become 
massless. The contribution of the poles to the integral ( I3.16P depends on the residues (I3.24p . 
As Z„ is exponentially suppressed for large momenta, the hole effectively decouples. The 
physical reason is that it is a collective excitation. Such excitations appear at length scales 
longer than the typical inter particle distance ~ T. Neglecting the continuum contribution, 
the spectral density can be approximated as 

Pp(Po) ^ vr ((70 - pj)5{po - a+) + (70 + P7)5(po + n+)) (3.26) 



tTij ^ aT : In this case medium corrections to the spectral density are negligible. Then 
no resummation is necessary and it is justified to use the free spectral function as a zeroth 
order approximation. It can easily be derived from (15.80 . 

p^-^(po) = 27r(2/+ m,)sign(po)5(/ - m^). (3.27) 
3.2.2 Computation of the Relaxation Rate 

The computation of the relaxation rate from couplings to fermions is rather involved. We 
restrict ourselves to the cases nXj ^ aT and m.j <^ aT previously discussed. 

xrii ^ aT : Insertion of ( I3.27p into ( I3.16p in the simplest case rtXi =m2 = m<m/2, 
|q| = yields 

/ /2m\2\^/^ 

ro = — m 1- — ) (l-2/^(m/2))^(m-2m), (3.28) 
87r \ \ m / J 

where we have used the identity f^^{uj)fF{po)fF{^ — Po) = 1 — /f(po) — fpi^ — Po)- The 
result ( I3.28P can be compared to the contribution (14. 2 p from decay into scalars. The 
kinematical restrictions are the same, but while the decay into scalars is enhanced by 
induced transitions at high temperature, for fermions in the final state there is a suppression 
due to Pauli blocking. 



tTii < aT : Insertion of iK26\) into flXTG]) yields for w > 0, |q| = 



Imn^(a;) = -—u' ( 1 " ( ^ ) ) (1 " 2/f(u;/2)) e{uj - 2071/), (3.29) 
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where we approximate 0+ ^ (p^ + dJtjy^^. This holds for |p| ^ 971/. This leads to 

To = f 1 - f ^ j j (1 - 2/i.(m/2)) ^(m - 2071/). (3.30) 

When is very heavy and the \E'j form a relativistic plasma, m ^ 9}t/ ^ rrij, the momenta 
of the decay products are large, |p| ~ m/2 ^ VJlj and f l3.30p can be used to compute 
the relaxation rate. Equations fl3.28p and fl3.30p show that in these limits the thermal 
masses qualitatively act as kinetic masses and close the phase space for the decay near 
Tc ~ m/(2a). Despite of this, the correct quantitative dependence of the relaxation rate 
on T even in this regime cannot be reproduced by replacing intrinsic masses with thermal 
masses, as one can easily see by comparing the exponents of the first brackets in (I3.28P 
and (ICTj) . 

Furthermore, the expression (I3.30p of course receives corrections near T^, where the 
involved quanta have momenta |p| < OJt/ and the details of the fermionic spectrum 
have to be considered. However, in the regime where they would be relevant the rate is 
suppressed due to the effect of Pauli blocking. The /i;'(m/2)-term already leads to an 
effective suppression at T ~ m <C m/(2a) for a ^ 1, thus it is likely that at Tc ~ m/{2/a) 
contributions from processes with bosonic final state typically dominate. This also implies 
that the closure of the phase space due to large thermal masses discussed for the trilinear 
scalar interaction is practically less important for the decay into fermions. The effect can 
be seen in figure [91 

Finally we would like to emphasise again that these conclusions are valid in the quasi- 
particle regime only. As discussed in section 13.1.31 and shown in figure [8] for the trilinear 
interaction, contributions from processes that are formally of higher order in the couplings 
can be of the same size as the leading order processes at high temperature. The suppres- 
sion due to Pauli blocking that is observed in the quasiparticle regime is less efficient in 
this regime because diagrams of higher order also contain processes with bosons in the 
final state. That can be seen easily by cutting the diagram shown in figure [10] through 
the scalar self-energy insertion. When there are bosons running in the loop, this cut yields 
contributions from scatterings with only bosons in the final state. 



3.3 Quartic Interaction 

The leading order contribution to Imll^ from the quartic interactions /ij0xf/4! comes from 
the setting sun diagrams shown in figure [T]l) and e). Analogue to f l3.2p it can be written 
as 

Imn;(c.) = -'^fs'i^) I (^^^p(^>o)A<k(fco)A<„p_,(a; -po - ^o)- (3.31) 

Again, the integrand is proportional to a product or spectral densities, in this case 
PipiPo) Pikiko) Pid-p-ki^ — Po ~ ko). In complete analogy to the discussion in section 13.11 
one can use finite temperature cutting rules to interpret the discontinuities in terms of 
elementary processes. As for the trilinear interaction, if the Xi resonances have a narrow 
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Figure 6: Processes contained in the setting sun diagram in figure [Hi), solid lines represent 
(f) propagators and dashed lines xi propagators. 



width and well-defined approximate dispersion relation, the integral f l3.3ip only receives a 
significant contribution if these processes are possible for real quasiparticles in the initial 
and final state. However, there is a significant difference as the setting sun diagrams at 
finite temperature contain scattering processes (pXi XiXii see figure |6l In contrast to 
decays, these processes are kinematically always possible, though at low temperatures they 
are suppressed by the statistical factors, reflecting the lack of scattering partners at low 
density. With (13.41) it is straightforward to derive the leading order expression 

11 J (ivrjy SWipWikWii 

X (^((l + /,)(l + /2)(l + /3)-/i/2/3) 
+ (/l(l + /2)(l + /3)-(l + /l)/2/3)) 

{5{uj + uJip - - uJi\) - 5(a; - cuip + cu^k + cuji)) 

+ ((l + /0/2(l + /3)-/l(l + /2)/3)) 

{5{u - ujip + Wik - Wji) - 5{u + Uip - ojik + c^ii)) 

+ (/l/2(l + /3)-(l + /l)(l + /2)/3)) 

{5{uj + Wip + Wjk - tUii) - 5{uj - ujip - Wjk + Wii)) ) , (3.32) 



here written in a symmetric way with the notation Wip = (p^ + m^)^/^, /i = fsi^ip) 
and so on. Again, the 5-functions ensure energy conservation and the statistical factors 
satisfy detailed balance. The first line in the brackets comes from decays and inverse 
decays (j) ^ XiXiXii which are kinematically forbidden for m < 3mj. The remaining lines 
correspond to scatterings ipXi ^ XiXi- Kinematically they are always possible, but as they 
require bath quanta in the initial state, they are suppressed by the occupation numbers 
appearing in the statistical factors and vanish in vacuum. At high temperature they provide 
the dominant channel of dissipation. Using resummed perturbation theory - replacing the 
internal lines of the diagram in figure [TJi) and e) by resummed propagators as in [T})) and 
[7b) - the dissipation rate at T mj can be approximated by [69] 



To-E^^- (3.33) 
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Figure 7: Feynman diagrams for the self interaction, the gray circles indicate 

resummed propagators. 



4 Comparison of Relaxation Mechanisms 
4.1 Decays vs Scatterings 

It is instructive to compare the contributions from scatterings and decays to estimate how 
efficient a heavy particle </> can dissipate or be produced from the plasma around T ^ T^. 
We focus on the trilinear and quartic scalar couplings because the decay into fermionic final 
states is suppressed by Pauli blocking. To determine Tc, specification of C^-i^t is required. 
For simplicity we chose £^.int = The temperature dependent dispersive part of 

the self energy to leading order comes from the diagram shown in figure [7^) 

= ml + A, / ^ ^2 hT\ (4.1) 

In the simplest case with m\ = m2 = m^, hi = h2 = h and Ai = A2 = A, (13. 6p for |q| = 
the rate (13. 6p reduces to 

' (l + 2/B(A//2))^(a;-2Mj, (4.2) 

where we have replaced the intrinsic mass by the thermal mass M^^, but neglected the 
thermal correction to the mass m because of the weak coupling. The scattering rate is 
computed from the diagrams in figure [T]i) and e) with resummed Xi propagators and can 
be found from I K33^ . At ^ (6(m2 - 4m2)/A)^/2 ^^kes the value 

, h'^m ( /2mv\2\ 
^•(^c)^2— — 1- — ^ . (4.3) 




scatter \^ Ci — ir>o \ i- \ 

1287rA V V m 



This can be compared to the decay rate (14. 2 p at T = 0, FdccaylT" = 0) and the maximal 
decay rate. The temperature Tmax at which (14. 2 p is maximal fulfils the equation 



(6((2mj2 - vr?) + T^^A) + Tl,Asinh(m/T^,.) = (4.4) 



^^The thermal mass correction here is ^ X^^^T rather than ~ AT because it arises from a local (tadpole) 
diagram, see figure [7^). 
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as one can easily see by demanding dT decay /dT = at T = Tmax- If Tmax ^ the sinh 
can be Taylor expanded and 

Tmax - ij{m^-i2m^f)\ • (4.5) 
With M^{Tmax) ~ {m^ + (m/2)2)/2 one can estimate 

raecay(T_) - ^a/I f 1 - (4.6) 



iTTm V A \ \ m 
where we have expanded fBinT-) as Tmux ^ fn- This allows to compare 

Tdecay (Tmax ) 
rdecay(^ = 0) 

■^scatter (-^c) 
rdecay(r = 0) 

r scatter (Tc) 



decay 



(T„ 




(4.7) 
(4.8) 
(4.9) 



Remarkably (14. 9 p does not depend on m^. In the interesting case <^ m the ratios 
are very simple. The parameter that governs the maximal amount of amplification of the 
decay rate by induced transitions is ~ (24/A)^''^. The parameter that determines whether 
relaxation via scatterings is efficient at Tc is (/i^/ X){m/g)'^. For g/m ~ h the rate rscatter(^c) 
can easily be bigger than rdccay(^ = 0) or even rdccay(^maa;) if A is sufficiently small, where 
one of course has to keep A ^ /i to be consistent with the initial assumption that the 
interactions that thermalise the bath are stronger than those that couple it to 0. 



4.2 The Importance of Off- Shell Contributions 

The importance of off-shell processes, which appear to violate energy conservation in the 
quasiparticle picture, can be parameterised by the widths Fj of the bath fields. At large 
temperatures this effect can be stronger than naively expected. 

A quantitative evaluation of fl3.15p requires specification of the interactions C^-i^t in 
( 13. ip that keep the bath in equilibrium. The dominant contribution to the thermal mass 
of scalars Xi can come from local (tadpole) diagrams which give purely real, momentum 
independent contributions to the self energy. For a AjXi/4! self-interaction these come 
from the diagram shown in figure Uh)- The resulting mass shift is given by (14. ip . Mj ~ 
(m?-|-AjT^/24)^/^ ~ A|^^T. In order to study the effect of dissipation via off-shell processes, 
we also have to include the width of the ^-resonances. The leading order contribution to 
the Xj-width from such interaction appears only at second order, see figure Wp), and is 
much smaller, cf. fl3.33p . Thus, the leading order contribution to the width may come 
from another coupling. 
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Figure 8: Dissipation rate Fq for the zero mode of from trilinear interaction computed 
with bare propagators - dotted blue, from equation fl3.6p - in comparison with the result 
obtained using resummed Xi propagators with thermal masses only (black dashed, from 
equation (13. 6p with masses replaced by thermal masses) and thermal masses as well as 
widths (red dotted, numerical evaluation of (13.151) ). In the upper row we chose a = 0.1 
(Fi/Mi ~ lO"'^), in the lower row a = 1 (Fj/Mj ~ 10 ^). The left column corresponds 
to «! = ^2 = 0.2 (Ml = M2 = O.IT), in the right column we chose a\ = 0.1, 02 = 0.3 
(Ml = 0.05T, M2 = 0.15T). ). Vacuum masses are small in all cases, mi = m2 = 0.1m. 



Therefore we in the following distinguish one coupling constant a that appears in the 
thermal masS Mi — OiiT/2 and another coupling cui that appears in the width. Equation 
(I3.15P shows that the inclusion of off-shell processes requires knowledge of the functional 
dependence of the discontinuities on energy and momentum. Thus simply taking F^ = afT 
is not sufficient. This also requires specification of the interaction that gives rise to the 
Xi-widths. We chose a Yukawa coupling to some massless fermion in the bath. That 
allows to use the formulae given in section IX^ with Xi as the external particle. This choice 
gives a good numerical stability when evaluating (I3.15p . Results for different parameter 
choices are shown in figure O 

^°For the interaction one can identify = (Ai/6)^/^, see ()4.ip . 
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Figure 9: The different rates in comparison: trilinear (black dashed), quartic (blue dotted) 
and Yukawa (red solid) ioi Y = h = gjim^pX) = 10~^ A = 1/24, a = 1/12. For this 
parameter choice the rates from Yukawa and trilinear interactions are equal at T=0. For 
the quartic interaction we only considered the contribution from scatterings ( I3.33p . 

5 Relaxation of a Fermion 

We now proceed to the relaxation of a fermionic field that is weakly coupled to a thermal 
bath. 



5.1 Fermionic Nonequilibrium Correlation Functions 

As in the scalar case there are two independent two point functions ^^^(xi, 0:2) and 



S>p{xi,X2) = (^a(a;i)^/3(x2))c (5.1) 

5<^(Xi,X2) = -{^l3ix2)^aiXi))c, (5.2) 

see also f lA.3l) and (IA.4I) . We again define a spectral function and statistical propagator 



S (Xi,X2) = i{S^{Xi,X2) - S^{Xi,X2)) 
S^{xi,X2) = ^ (5'>(xi,X2) + S'<(xi,X2)) 



(5.3) 
(5.4) 



Their time evolution is governed by the Kadanoff-Baym equations (1A.38P and ( ]A.39p . The 
solutions for Majorana fermions have been found in P5|I24]. With the symmetry relations 
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given in appendix lA. 1 1 the generahsation to Dirac fermions is straightforward'^^!. 

S"AU-t,)=i -e-(*-*^)p,(a;) (5.5) 



oo 2vr 



and 



+ f ' dt'S^ih -t') r dt"J:+{t' - t")S^{t" - t2) (5.6) 
Jo Jo 



with 



As in the scalar case, the spectral density for the nonequilibrium field coincides with 
the equilibrium expression (13.171) even for large deviations from equilibrium as long as the 
bath is sufficiently large. Fermion spectral functions are, as pointed out previously, in 
general rather complicated. Usually replacing vacuum masses by thermal masses does not 
even in the quasiparticle regime give a consistent approximation to the propagators. We 
are interested in the behaviour of a massive fermion that is weakly coupled to a thermal 
bath, thus it is consistent to consider the situation m ^ YT only. Then the effect of the 
medium on the dispersion relation of is small and can be neglected with respect to the 
vacuum mass m. The thermal correction to the width is, however, not negligible as there 
is no large width in vacuum in comparison to which it could be neglected. We assume that 
the discontinuity of the fermion self energy is a pure vector, Sp(w) = 2i{aci{u);^+ bq{ijj)i{) . 
Here u is the four velocity of the bath, which in its rest frame is simply (1,0) and aq(ci;) 
and bf^{u) are functions of g = (w, q) that depend on the interactions of This of course 
excludes a number of interesting applications such as CP violating terms, but the inclusion 
of those is straightforward and does not affect the kinematical aspects we are interested in 
here. Then (I3.17P can be approximated by a Breit Wigner function, 

p^iuj) ^ 2(g/+ m)^^ \ (5.8) 

where cUp = (m + q^)^''^ and 

Tq = 2 (a^— + . (5.9) 
The u integration in (15. 5p can, analogue to (I2.20p . be evaluated using Cauchy's theorem, 
S-ih - t,) ^ e-^-^"^-'^'/^ (^^70 cos(a)q(ti - t,)) - ^^^-^ sm{uo^{t, - ta))) • (5.10) 

•^^I would like to thank Wilfried Buchmiiller, Alexey Anisimov and Sebastian Mendizabal at this point, 
in collaboration with whom the propagators for Majorana fermions were first found. 
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Figure 10: Fermion self energy, the lines represent \1/ (solid with arrow) and the fermion 
(dashed with arrow) and scalar (solid) in the bath. 

Fourier transformation and another use of Cauchy's theorem allow to evaluate the memory 
integral in the second line of ( I5.6P 

ZCUq 

X ((m -7q) cos(a)q(ti - ^2)) - «7o<^qSin(a)q(ti - ta))) (5.11) 

Together fIS.lOp and (15. lip provide approximate analytic expressions for the correlation 
functions (15. 5p and (15. 6p . We re-emphasise that it is in general not a valid procedure 



to obtain non-equilibrium propagators from the equilibrium correlation functions (1A.19 
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by replacing the Fermi-Dirac distribution function fp with some generalised distributio 
(15. 5p and (15. 6 P have been found as solutions to the full non-equilibrium equations of motion 
( ]A.38p . ( lA.39p . Analogue in the scalar case, all correlations are damped with a rate Fq, 



which can be interpreted as relaxation rate in accordance with [5l]- This interpretation 
holds even for large deviations of \l/ from equilibrium. 

5.2 Relaxation via a Yukawa Coupling 

When the dominant contribution to Fq comes from a Yukawa coupling of strength Y to 



a Dirac fermion \l/ and a real scalar (p in the batlo, it can at leading order be computed 
from the Feynman diagram shown in figure [TOl 

^^{uJ) = -^Y'f^\co) J (^5<(po)A>_^(po-c.) 
f d^p 

= -iy^fp^^) J ^^^/f(Po)/b(w -Po)pp(Po)Pq-p(w -Po)- (5.12) 



^^As in the scalar case on can of course define f{q,uj,t) = —S^{t;uj)p^^{uj) + 1/2, but this requires a 
dependence of / on time and ofF-sheh energies. 

■^•^This means we consider a Lagrangian very similar to p.ip : we could identify the out of equilibrium 
fermion with and the bath fields 4" and cf) with 'if 2 and cj), respectively. We prefer to use an alternative 
notation here in order to avoid confusion as (jj in the previous sections has always been a weakly coupled 
field out of equilibrium while (f> here is a bath field and 4* is out of equilibrium. 
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Here 5*^ and are the thermal Wightman functions of the fermion and scalar in the bath, 
see ( lA.igp and (lA.lSp . As in the scalar case, they have to be resummed to obtain consistent 



results. The relaxation rate is again given by a product of spectral densities. Once again, 
the spectra in the bath may be highly complicated, but as long as the contributions from 
narrow poles dominate the integral fl5.12p . the discussion from section [STT] can be repeated. 
The poles only contribute if quasiparticle decays (case (a)) or Landau damping (case (b)) 
are kinematically allowed on-shell. A Yukawa coupling is a three body interaction, thus a 
kinematically forbidden region (c) in which none of the involved quasiparticle can decay 
into the others is expected. In this region Fq only receives contributions from processes 
involving intermediate off-shell quasiparticles and is suppressed. The phase space closure 
is more relevant in this case than for the decay of scalars into fermions because Pauli 
blocking does not suppress the efficiency of the decays near the critical temperature as one 
of the final state quasiparticles is a boson and the /^-term in fl5.12p can compensate the 
suppression by the fp-term. 

In section 13.21 the two cases m ^ aT and m <C aT were considered, where a is a 
typical coupling for the fermion in the bath. Here we focus on the case m ^ aT. Then 
no resummation is necessary for the fermion propagator. For scalars replacing the vac- 
uum masses by their thermal counterpart is a consistent approximation for the resummed 
spectral density if the only relevant resonances are the screened one particle states, their 
widths are small and the dispersion relation is approximately parabolic. This is e.g. the 
case if the leading order contribution to the effective mass comes from tadpole diagrams 
as shown in figure [TH) or a Yukawa coupling (67] . Then we can use the expressions 

Pq_p(a;-po) ~ 2TTsign{uj - po)6{{q - - M^) (5.13) 
Pp(po) ~ 27rsign(po)(y+m)5(p2-m2) (5.14) 

to evaluate (I5.12p . where the scalar spectral density p is resummed, M contains vacuum 
as well as thermal contributions while m is a chiral vacuum mass. After performing the po 
integration reads (cf |51) ) 



S-(a;) = -^r2 



(27r)2 2lii2fi2 

(^{1 - fF{Q,) + fB{n2)) ((1^170 -P7 + ni)6{u -Q,-Q2 

+ (MQ,) + /b(1^2)) ((^^i7o -P1 + m)5(a; - + Q^) 

+ {Qi^o +P7 - + Q1-Q2)] ] (5.15) 



with Cli = (p2 + rh^)^ and = ((q — p)^ + M^)^. The kinematic 5-functions in (I5.15P can 
be compared to those in (13. 5p and lead to similar kinematic restrictions and regions (a), (b) 
and (c) in which exchanges energy with the bath via on-shell decays and inverse decays 
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Figure 11: Vq in units of its zero temperature value as a function of T/m for |q| = m/5. 
The analytic expression obtained in appendix lA. 21 without (blue dotted) and with thermal 
mass (red solid) for the scalar is compared to the approximation ( I5.16P with thermal mass 
(black dashed). For the scalar we took a vacuum mass rh = 0.4m. and a quartic self coupling 
of strength A = 1/5. 

in case (a), on-shell Landau damping in case (b) or only off-shell processes in case (c). 
As in (13 ■5p . the integral in ( I5.15P is exactly zero in case (c), but would receive a non-zero 
contribution if the widths of the resonances in the bath were taken into account. 

When thermal corrections to Pp are relevant, also the fermion propagator has to be re- 
summed. The evaluation of Fq then becomes more involved, but it is obvious that, when the 
widths are small, (I5.12p is again dominated by pole contributions and kinematic arguments 
based on approximate energy conservation in reactions between on-shell quasiparticles can 
be used. However, as usual for fermions, this in general cannot be done by simply putting 
in thermal masses by hand. The decay into fermions with resummed spectral densities 
(I3.18P has e.g. been studied in [5l|6^ where it was found that Fq is non-zero above the 
critical temperature associated with the decay into particles because phase space for the 
decay into scalars and on-shell holes closes at higher T only. 

In the following we will simply assume that ( I5.13p . ( I5.14p hold at all temperatures 
of consideration. If m is negligible the integral in fl5.15p can be solved analytically for 
arbitrary q and uj, see appendix IA.2[ The analytic expression is rather complicated, but 
for nonrelativistic on-shell quasiparticles with |q| < 0.5m the quasiparticle width Fq can 
be well-approximated by 



with m-t = (m^ ± Af^)/m. For |q| ~ m this approximation becomes inaccurate. Figure 
[TT] compares the approximation (I5.16P to the analytic solution of (I5.15P given in appendix 




(5.16) 
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[Q The forbidden region (c) here shrinks to a point because one of the particles in the 
loop is massless, thus qthi = gth2- 

Note that the requirement to keep the fermion in the loop massless is indeed quite 
strong. Even if the fermion is massless in the vacuum, m = 0, it will in general obtain a 
thermal mass via radiative corrections because thermal masses are not protected by the 
chiral symmetry |28| . see fl3.18p . 

6 Comparison to Boltzmann Equations 

Many nonequilibrium phenomena, in particular in cosmology, can be studied in terms of 
Boltzmann equations. They are by construction Markovian and formulated in term of 
number densities or single particle phase space distributions for classical particles. There 
are phenomena, including coherent oscillations or memory effects, that cannot be described 
in this framework. In absence of these one should recover the Boltzmann equations as a 
limit of the quantum field theoretical description. 

The rate Fq determines at which rate the field mode q gains or loses energy. In the 
quasiparticle regime, the gain and loss rates can be understood in terms of approximately 
energy conserving decays and scatterings of quasiparticles. This suggests that the system 
can be viewed as a gas of quasiparticles and one is tempted to simply replace m — )■ M or 
Wq — 7- fiq in the standard Boltzmann equations. In the following we show that this is not 
necessarily correct. 

We consider the simplest case, a single real scalar field coupled to a large thermal bath. 
We restrict ourselves to the quasiparticle regime as it is clear that outside this regime Boltz- 
mann type equations have no meaning. For a single scalar without additional quantum 
numbers (such as flavour), coherent oscillations play no role. When the bath tempera- 
ture is constant, the quasiparticle approximation allows to sum up all memory effects and 
obtain the analytic expressions fl2.20p and f l2.2ip for the spectral function and statisti- 
cal propagator. We also assume that the dispersion relations are in good approximation 
parabolic. 

The Boltzmann equations are formulated in terms of one-particle phase space distribu- 
tion functions. As the definition of a number density in the field theoretical approach is 
ambiguous, these are not well-suited for a comparison between the two. We will therefore 
perform the comparison in terms of the energy stored in each mode, which is a well-defined 
observable in both approaches. 

Boltzmann Equations Let us consider the Boltzmann equation for particles of mo- 
mentum q and energy iJq. The competition between a gain and a loss term, 7q and 7q 
respectively, determines the change of the particle number /q(t) for the mode q, 

dtUit) = {1 + Ut)h< - Uth> . (6.1) 
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When the medium is in equihbrium, production and decay rates satisfy the detailed balance 
relation, 

7< = e-^^n> = /B(i?q)7q , (6.2) 

which implies 7q = 7q — 7q , analogue to ( I2.24p . The rate 7q is usually computed at T = 
from scattering cross sections or, via the optical theorem, the imaginary part of the self 
energy, 

Imn^(Eq) 



7q 



T=0 

Using these relations, the Boltzmann equation (16. ip can be written in the forr^'^^' 



(6.3) 



dtUt) = -7q(/q(t) - fB{E^)) , (6.4) 

with the solution 

Ut) = UE^) + (/q(0) - UE^)) e--^-' . (6.5) 
The energy for each mode is obtained by multiplying (16.51) with E^ 

iMt) = En [fsm + (/q(0) - UE^)) e-^'^') . (6.6) 

Kadanoff-Baym Equations In the field theoretical description the energy per 0-mode 
is given by (12.70 . 



■■t2=f 



The (0) (0)-part corresponds to the energy that is stored in the mean value (0) of the field 
while the A+-part comes from its fluctuations, to be interpreted as (quasi)particles. An 
explicit expression for e^(t) can be found by inserting (I2.2ip and (I2.22p . 

<(^) - if Ksin^(fi.^) + fiScos^(fiqt))e-^<^* 



21]2 

q 



^q,mY^q,m 



Sin(fiqt) COs(f2qt) — f2q) 6 



K2 



A+. 

q,m 






2 






A+. 

q,m 




"q 


2fi2 


\ 2 



' ""q ^ g-Fqt 



^*Note that the argument of fq{t) is a time while fsiEq) does not depend on time and we have exphcitly 
written the particle energy as the argument. 
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The first three terms represent the energy of the mean field (0). We will ignore them in 
the following and set ^q^n = 0q,in = 0. 

Comparison Obviously neither identifying with Uq nor with generally leads to 
equivalence between (16. 6p and (16. 7p . The oscillations in (16. 7p . which contains the solution 
to a second order differential equation, remain present in the limit f2q — cjq, Fq — )■ 7q. 
Therefore the system cannot even be described by Boltzmann equations if the relaxation 
rate Fq is identified with its vacuum value and all modifications of the dispersion relation 
are neglected. 

However, a consistent comparison between a quantum mechanical and a classical observ- 
able can only be done if the quantum system is set up in an initial state has a counterpart 
in the classical theory, namely one of definite initial (quasi)particle number. 

We first consider the simplest case and assume that initially is in equilibrium, but at 
a different temperature T than the bath. One can then think of the system as a two heat 
baths at different temperatures T and T that are brought in touch at t = 0. One bath (0 
only) is much smaller than the other, neglecting backreaction the whole system is expected 
to equilibrate at temperature T for times t 2> 1/Fq. The initial density matrix can be 
written as Qmi oc e~^^'^ ® e"''^''^'*'' and the initial statistical propagator for is given by 

A+(ti, t2) = ^ Q + fBi^Sj cos((ti - t2)nq) at U (6.8) 

Here fs is the Bose-Einstein distribution with temperature T. The energy (16. 7p then reads 

•)2 I , ,2 



e 



i + MnSj + (/B(f^q) - /B(l^q)) e-^''*^ . (6.9) 



2% 

The solution (16. 9p describes the transition from one thermal state to another. It can be 
compared to (16. 6 P with /q(0) = /^(i^q). Apart from the additional term 1/2 in (16. 9p . 
which is simply the vacuum energy of the mode q, the expression in the bracket looks 
exactly like the solution for a Boltzmann equation for quasiparticles with energy f2q and 
damping rate Fq. However, it is not multiplied by Q^, but {Q'^ + uj'^)/{2Qq). If one loosely 
interprets the term in the brackets as quasiparticle number for mode q, it appears that it 
follows a Boltzmann equation while the energy does not. 

The same effect can be observed when inserting the more general initial conditions 

q,m n V2^^q'ini;) ^q,m q,m ^ \2 ^ J 'i'^^^ J ^ 
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sometimes referred to as tsunami initial conditions. Here /q,ini is a general function of q 
that can loosely be interpreted as initial occupation number of that mode. The solution 
in this case reads 

4{t) = (^Q + Mn^)j + (/q,M - fB{%))e-^^'j ■ (6.10) 

Again, the term in the brackets may be interpreted as the number of quasiparticles per 
mode q that follows a Boltzmann equation. However, the energy, which is the observable 
quantity, does not. Obviously (16.91) is a special case of f lG.lOp with /q,ini = /^(fiq)- One 
way to interpret this result can be seen by rewriting it as 

ej(t) = fiq {fsin^) + (/q,ini - fsi^)) C^^'^') 

The first line is a Boltzmann equation for quasiparticles. The second line may be inter- 
preted as a vacuum energy. However, it cannot be removed by simply shifting the energy 
of the ground state by a constant because it depends on temperature and time. It appears 
even tough we have not considered the coupling of the particles to the time dependent 
background (0) , which is of higher order in the coupling between (j) and the bath. There- 
fore, even though the exchange of energy between different field modes and fields can be 
understood in terms of elementary processes involving quasiparticles, our result suggests 
that the total energy of is not that of a quasiparticle gas. A similar observation was 
made in equilibrium in [6], where it was found that due to this effects the contribution of 
to the total pressure of the system can be negative. 



7 Application to Reheating after Cosmic Inflation 

We now discuss an application of the previous results. It has been suggested in j32] 
that large thermal masses in the primordial plasma can impose an upper bound on the 
temperature in the early universe. The authors argue that, if reheating is fuelled by the 
decay of infiaton particles, this process would stop when the thermal masses of the would-be 
decay products become sufficiently large to make the decay kinematically impossible. 

Any bound on the temperature in the early universe has far reaching cosmological 
consequences as it determines abundance of thermal relics. It is, for instance, of crucial 
importance in the context of the gravitino problem [23] and thermal leptogenesis |34| . A 
large temperature can also be related to the fate of moduli |35) and the possible decom- 
paticfication of extra dimensions |36] . 

7.1 Cosmic Inflation and Reheating 

Cosmic inflation jTD], the idea that the universe underwent a period of accelerated expan- 
sion during its very early history, has proven an extremely successful concept in explaining 
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various unsolved mysteries of big bang cosmology. These include the flatness and horizon 
problems, the absence of topological defects and the origin of the temperature fluctuations 
in the cosmic microwave background, thought to be the seeds of structure formation pGj. 
Accelerated expansion requires a negative pressure component to dominate the energy con- 
tent of the universe. The simplest way to realise this is by the potential energy V{(j)) of a 
scalar field 0, the inflaton. During the inflationary era all other forms of energy are diluted 
and virtually all energy is stored in the zero mode of the mean fleld (0) , which slowly rolls 
down the potential V{(f)) until it starts oscillating around its minimum. It leaves a cold and 
empty universe with very small entropy. The conditions at the beginning of the radiation 
dominated era, corresponding to a hot big bang, are then created by reheating. During 
this process (0) performs damped oscillations around its potential minimum, dissipating 
its energy into all other degrees of freedom. This flUs the universe with particles and heats 
up the primordial plasma. 

We restrict ourselves to scenarios where only one scalar inflaton field (j) is relevant during 
reheating. This includes scenarios in which the underlying physics can be parameterised 
in terms of a scalar degree of freedom. We describe the interactions of (p with matter by 
the Lagrangian fl3.ip . where (p is identified with the infiaton. Xi are some scalar 

and fermionic fields that directly couple to the infiaton and via C^^i^t and £.i,^int directly 
or indirectly interact with all other fields in the primordial plasma (including the SM). 
This can be regarded as a toy model, but is sufficient to describe the kinematic aspects of 
dissipation into bosons and fermions we are interested inl^. 

As soon as (j) has dissipated some fraction of its energy into other degrees of freedom, 
the presence of the plasma formed by its decay products affects the ongoing reheating 
process, which in consequence has to be treated by methods of nonequilibrium fleld theory. 

The details of the reheating process are largely unknown due to our lack of knowledge 
about the underlying physics. It is nevertheless possible to study general features and 
their cosmological implications. Here we study the effect of thermal masses in the plasma 
on the time evolution of the temperature. We are in particular interested in two quanti- 
ties, the maximal temperature Tmax that the universe was exposed to and the reheating 
temperature Tr at onset of the radiation dominated era. 

7.1.1 Effective Masses 

While the part of the inflaton potential that is relevant during inflation can be probed 
by cosmological data |71 | . very little is known about the region around the minimum 
and the inflaton mass iWlfj during reheating. Even without thermal effects M should be 
regarded as an effective mass that may depend on the values of other flelds. We make 
two assumptions about M. First, it is much larger than the vacuum masses rrii, nXj of the 
other flelds. Second, it does not depend on time during reheating. This is consistent with 

^^It is likely that gauge fields have played an important role during reheating, but the kinematic aspects 
we are interested in are the same as for scalars. 

^^Here we use a capital letter to denote the mass of at T = because it should be regarded as an 
effective mass that might depend on the values of other fields. 
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Figure 12: Contribution to the Xi self energy from coupling to (0). 



the assumptions that is the only relevant out-of-equilibrium field and weakly coupled, 
making thermal corrections to its dispersion relation negligible. Then we can treat the 
0-mass as a constant parameter M. 

Regarding the effective masses and dispersion relations of the other fields, there are 
three different contributions. The first is their vacuum mass. The second comes from 
interactions with the plasma. If the plasma is in equilibrium it can be treated by the 
methods discussed in the previous sections. The third is due to the coupling to the mean 
field {4>). For the Xi if arises from diagrams as shown in figure [T2l 

In the following we always neglect all vacuum masses except the </)-mass M. This is 
justified since Tr is likely to be far above the electroweak scalqfj. We furthermore assume 
that all medium effects on dispersion relations in the plasma can be parameterised by 
momentum independent thermal masses. This assumption is for simplicity and due to 
our lack of knowledge about the exact decomposition and interactions of the primordial 
plasma. It can, however, also be justified when 0-quanta decay into much lighter states 
with large momenta (cf. sections 13.21 and 14. ip . Following the discussion in section ISTT] a 
consistent treatment requires to use full resummed propagators in the bath. This includes 
not only the thermal masses, but also the widths. 

For a quantitative treatment of both, the interaction terms C-^.-mt and ^iii^int have to 
be specified. Following the considerations in section 14.21 we again introduce two different 
couplings a and a that govern the real and imaginary part of the self energies. The 
thermal masses can be estimated as Mx — aT/2. Equation fl3.15p shows that a consistent 
treatment of contributions to the relaxation rate involving virtual Xj, quanta requires 
specification of the functional dependence of the discontinuities on energy and momentum, 
thus taking Fj ~ a^T is not sufficient. As in section 14. 2[ we again assume that the main 
contributions to the Fj comes from Yukawa couplings to some fermions in the bath. This 
allows to use the result from section 13.21 

Finally, there is the contribution to the effective masses from diagrams as shown in 
figure [T^ These terms mix Xi and X2- For mi,m2 <^ aT the effective Xj-mass after 
diagonalisation of the mass term is 



Here we have for simplicity only considered the trilinear interaction for the Xi- A similar 
term ~ ^(0) arises as correction to the fermion mass. 

^^It is clear from the previous sections that the inclusion of vacuum masses is straightforward. 



r^a^T^ + g{(P). 
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7.1.2 Applicability of our Results 



Denoting by $ the amplitude of (0), (17. ip allows to distinguish the regimes > a^T'^ 
and < a^T^. The expectation value (0) = (0q(t)) oscillates rapidly at frequencies 
> M, faster than the relaxation time in the plasma. Therefore the formulae used in the 
previous sections, valid for coupling to a large bath in equilibrium, are not applicable 
for (7$ > a^T^. In this case all degrees of freedom are far from equilibrium!^ and there 
are various nontrivial collective effects such as parametric resonance, which can not only 
produce particles with masses larger than M, but also potentially observable relics [72 | [73 ] . 
However, we are interested in the effect of thermal masses (contributions from coupling to 
the background plasma). It is clear that in this regime they only give rise to a subleading 
correction. 

For (7$ ^ a^T^ the only relevant timescales are the relaxation times of the decay 
products (tx) and (p (r^). The relaxation times tx are usually much longer than the 
period of one </)-oscillation, but short compared to t^, which characterises the duration of 
the reheating process. On timescales < the system can effectively be described as a 
single scalar field that is weakly coupled to a thermal bath with many degrees of freedom 
and the methods from the previous sections can be employed. In this regime dissipates 
via perturbative decay of single quanta. 

Since T and $ change with time, the system can pass through various stages. In many 
scenarios the infiaton releases much of its energy during an initial phase of preheating 
during which the time dependent mean field dominates. However, the decomposition of 
the primordial plasma at the onset of the radiation dominated phase may nevertheless be 
determined by the temperature during a later phase of reheating by perturbative decay. 
The reason is that the density of any particles produced during preheating is diluted by 
the expansion of the universe during reheating. We denote by the energy density of 
while ept refers to the energy density of all other (7* relativistic degrees of freedom. We refer 
to the period during which > e^j as reheating stage. For < ept the universe enters the 
radiation dominated era. The temperature Tr at = is the reheating temperature. 
It is generally lower than the highest temperature Tmax that the universe was exposed 
to, see figures [MllTSl Indeed, during most of the reheating stage the temperature actually 
decreases while en/e^ grows. Since the 0-oscillations are fast, (averaged over a few 
oscillations) redshifts like nonrelativistic matter even at temperatures T > M. In contrast 
to that the bath quanta are relativistic and redshift like radiation as their average momenta 
in the plasma have values ~ T while thermal masses are ~ aT, leading to a relative dilution 
of particles produced at early stages. This makes Tr a very decisive quantity. 

We restrict ourselves to initial conditions corresponding to perturbative decay, including 
the possibility that this follows an initial stage during which other mechanisms dominated. 
As $ decreases with time the inequality (7$ > a^T^ remains valid at all later times. The 
question whether large thermal masses can stop the heating of the plasma by blocking the 
phase space then arises once aT ~ M. 

^^This also means that there is no well-defined temperature T in (|7.ip . but there is still a contribution 
due to the coupling to the plasma at order a^. 
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Figure 13: Example for a contribution from coupling to the mean field to the self energy 
n^, dashed and dotted lines denote Xi and X2 propagators, the crosses indicate (0). 

The relevance of the perturbative decay depends on the amount of energy still stored 
in when it takes over. We demand ^ while <^ a^T^. One can estimate 
^(p ~ ^(0) ~ M^$^ and epi ^ iT^g^:T^/30. This estimate for can also be obtained from 
(EZD with fiq = and A+.„ = A+,„ = A+.„ = 0, 

<W = ^(0lin + ^q0q,in)e-^'^*, (7.2) 

which for 0q <^ 0^ (slow roll) and q = reproduces the above estimate and is justified 
because is very weakly coupled and all energy initially is stored in the mean fielcj^. 
The usual estimate for is derived from classical thermodynamics using single particle 
distribution functions, see e.g. |Tj. Our results from section |6] suggest that it should receive 
corrections even in the quasiparticle regime, but we shall ignore this subtlety for the mo- 
ment. It straightforward to formulate the condition a/g ^ (7ry^)/(30M). For g^, ~ 100 
this means a ^ g/M which is generally the case. 

Equation (16. 7p . or f l7.2p . shows that the dissipation rates for the oscillating field (0) 
and its fluctuations near the potential minimum are both given by F = Fq=o computed 
in section I^FI . The frequency of the oscillations is r2q=o — M. It is not fully consistent 
to use fl6.7p or fl7.2p because the solutions for the propagators A~ fr2.6p . A"*" 02. 9p and 
field value (0) (I2.13p from which they were derived were found under the assumption 
of time translation invariant self energies computed from couplings to a large bath of 
constant temperature. This condition is violated in three ways during reheating. On one 
hand, Hubble expansion cools the plasma continuously. On the other hand, in the case of 
reheating backreaction (the effect that the dissipation of has on the bath) is certainly not 
negligible by the very definition of the process. Finally, contributions to the 0-self energy 
that involve itself (see figure [T3|) can have a sizeable effect despite the large number 
of degrees of freedom if (0) is sufficiently large. This last point in principle imposes the 
strongest restrictions on the applicability of our results, but is not relevant here as we 
restricted ourselves to initial conditions g^mi < a^Tf^^. 

Hubble expansion and backreaction can be parameterised by a time dependent tem- 
perature. As long as the relaxation time of the Xi is much shorter than the timescales 

^^Note that strictly speaking e^(t) is an energy per mode while is an energy density. 

'*°Equation (|2.22p . and thus (|7.2p . strictly holds only in the harmonic approximation of the potential 
y(0). In the late phase of reheating which determines this can be justified. Furthermore, the interpre- 
tation of Fq as damping can be extended beyond the linear regime |46) . 
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associated with ^-relaxation and Hubble expansion, there is a separation of timescales. In 
any case we are not interested in the details of the field oscillations. That allows to use the 
expressions for F obtained in section [2] with a time dependent temperature and interpret 
them as the dissipation rate at each moment in time. 

7.2 Thermal History during Reheating 
7.2.1 Boltzmann Equations 

In the following we use the results from the previous sections to study the time evolution of 
the temperature during reheating. We describe the system in terms of effective Boltzmann 
equations for (averaged over a few oscillations) and ej^. 

e^ + 3He^ + Te^ = (7.3) 
eR + 4HeR-re^ = 0. (7.4) 

Here H is the Hubble parameter. These equations can be written more conveniently in the 
variables $ = e^M~^a^ and R = eaa'^, 



with 

T = —[^r\ (7.8) 

where x = aM, a is the scale factor, Mp the Planck mass and the number of relativistic 
degrees of freedom in the bath. 



7.2.2 Dissipation into Bosons 

We first study the dissipation into the fields Xi- The relevant parameters are the inflaton 
mass M and the couplings g, hi = h2 = h. The fields Xi ^^d X2 can have different thermal 
masses if their dominant couplings ai and a2 to the rest of the bath are different. While 
keeping ai and a2 independent, we assume that their widths are governed by the same 
parameter ai = 0.2 = a. It would be straightforward to include the possibility of different 
widths, but this does not lead to any new effects (while Mi 7^ M2 can give rise to Landau 
damping as a new dissipation mechanism). Without loss of generality we take 02 > ai- 

There are two critical temperatures T^. = 2M / [ai+a2) and = 2M/\ai—a2 \ associated 
with the trilinear coupling, defined by the conditions M = Mi + M2 and M = \M2 — Mi\. 

*^This of course assumes that memory and coherence effects are not essential 
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The initial temperature Tj„j, which determines the initial value for e/j = vr^^f^T'^/SO, shall 
be smaller than T^, Tj„j ^ ~ M/a. We furthermore demand that is large enough 
that the energy stored in can potentially heat the plasma to temperatures higher than 



If (17.1 op is not fulfilled, thermal masses play no essential role because the effective Xr 
masses are dominated by coupling to the oscillating field (0). If (17. 9p is not fulfilled, the 
plasma never reaches a temperature at which they would be relevant. 

Consistent with all the above conditions we chose Tj„j = T^. = M/a, = 
Q.ha^TfJg = M^/{8g). We chose a = 0.1, hi = h2 = h, = 200 and M = WGeV. Then 
the only parameters to vary are the ratio g/M (which governs F/M as F oc g'^/M) and 
a (which governs T^/MyJ. These choices are for illustrative purpose and are not meant 
to resemble a particular inflationary model, the study of which we leave for later work. 
We are mainly interested in two quantities: The maximal temperature Tmax achieved and 
the reheating temperature Tr at the beginning of the radiation dominated era (defined by 
e_R = e<^)- We consider various representative choices for the remaining parameters. In all 
cases the dissipation rate for the trilinear interaction is computed by numerical evaluation 
of (I3.15P in each step in order to include off-shell corrections. For scatterings by the quar- 
tic interaction, which are always possible on-shell, we use the resummed rate (I3.33P and 
neglect off-shell corrections. 

I) /z, y -C g/M = 10~^, ai = ^2 = 0.2 : The results for this case are displayed in 
figure [m The quartic interaction is negligible. Thermal masses close the phase space 
for the decay — >■ X1X2 at Tc ~ 5M. Due to the equal thermal masses of xi ^ind X2, 
Ml = M2 = O.IT, Landau damping via processes Xi'P ^ Xj is kinematically never possible, 
see figure |H1 Off-shell effects are negligible for a = 0.1 and a = 1, corresponding to 
Ti/Mi ~ 8 ■ 10~^ and Ti/Mi ^ 8 ■ 10"^. Without resumming the Xi propagators (i.e. 
neglecting thermal masses) one clearly overestimates Tmax- However, for this choice of 
parameters, one would still get the correct value for Tr. 

II) h,Y g/M = 10~^, tti = ^2 = 0.2 : The results for this case are displayed in figure 
[T5l The situation is similar to case I), but T/H is larger and the plasma heated more 
efficiently. The widths Fj have negligible effect for a = 0.1 {Ti/Mi ^ 8 ■ 10"^), but allow 
a temperature Tmax slightly larger than for a = 1 {Ti/Mi ~ 8 • 10"^). Though it is 
conceptually interesting to see that processes involving virtual Xj-quasiparticles can sustain 

'^^This then also covers the case Tc < Tint < Tc because the plasma is eventually heated up to Tc- The 
case Tini > Tc is not interesting because Landau damping obviously allows heating to high temperatures. 




(7.9) 



For the effective masses to be dominated by thermal contributions we demand 



g<^im < a^Tli 



(7.10) 
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Figure 14: Temperature as a function of the scale factor for parameter choice I), M = 
lO^GeV, g/M = \/2- 10^^ and ai = a2 = 0.2. The plot in the lower right corner compares 
the evolution of T for a = 0.1 (Ti/Mi ~ 8 ■ lO""^, blue solid line), a = l {Ti/Mi ^ 8 ■ lO'^, 
red dashed line) and without resummation of the Xi propagators [Ti = 0, Mi = rrii <^ M, 
green dotted line). The other three plots show the temperature evolution (red line) and 
the relative contribution of (dotted line) and en (dashed line) for each single case: upper 
left a = 0.1, upper right a = 1, lower left without resummation. 



reheating for T > T^, this small correction is phenomenologically not relevant (especially 
given the uncertainty about the correct model for inflation). However, for g/M = 10~^ the 
beginning of the radiation dominated era falls into the epoch during which T is "frozen" at 
Tc. Therefore the thermal masses affect not only Tmax, but also Tr and the temperature 
during the radiation dominated era. 

Ill) h,Y <t:g/M = V2- 10-^ ai = 0.01, = 0.39 : The case that Mi and Ms are very 
different and a = 1 is interesting. Figure [16] shows that correct values for Tmax o-nd Tr 
can only be obtained using full resummed propagators, including thermal masses as well 
as widths. The reason is that for 0^2 3> ai the region (c) shrinks as Tf. — T^. — 4Mai/a|, 
see figure [T71 Then the dissipation by off-shell processes that only led to an insignificant 
increase of Tmax slightly beyond T^ in figure [15] can be sufficient to heat the plasma up to 
Tc. After it crossed the forbidden region, reheating can continue via Landau damping. 
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Figure 15: Same as figure [Ml but for parameter clioice II), g/M = \/2 ■ 10 



IV) g/M = h/2 = 10"^, tti = 0.01, aa = 0.39, a = 0.1 In this scenario we fix a = 0.1 
such that the Xi are in the quasiparticle regime (Fj/Mj ~ 8 ■ 10~^) and off-shell contribu- 
tions are negligible, but add the quartic interaction. The trilinear interaction cannot heat 
the universe to temperatures larger than T^. The quartic interaction mediates (pXi ^ XiXi 
scatterings at leading order. These are inefficient at low temperatures due to the low 
occupation numbers, but can become the main source of dissipation for large T. Fig- 
ure [18] demonstrates that scatterings mediated by the quartic interaction can allow for a 
temperature beyond Tc. 

7.2.3 Decay into Fermions 

The contribution from processes with fermions in the final state is generally thought to be 
smalQ While the dissipation rate into bosonic final states is amplified in a plasma due to 
induced transitions, for fermions there is a suppression due to the Pauli principle. These 
effects become relevant when /^ p'(M/2) is of order one, which is the case for T ~ M. 
The effects of thermal masses kick in around T ~ M/a. Thus, for a ^ 1 thermal fermion 
masses will not affect Fq significantly. In the quasiparticle regime this allows to draw 
two conclusions. First, perturbative decay into fermions is negligible at temperatures near 
the inflaton mass if there exists a coupling to bosons of comparable strength (F,^^^^ ~ 
F(^_^^^ at T = 0) due to the relative suppression F^^^i^/F^^^.^^ ~ (1 — 2fp{M/2))/{l + 

^■^It has been pointed out though that fermions may play a crucial role during preheating |74| . 
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Figure 16: Same as figure [TU but for parameter clioice III) witli g/M = a/2 ■ 10 ^ and 
ai = 0.01, as = 0.39. 



2/s(M/2)) = (tanli(M/3/4))^ ^ M^(4T)~^. Second, tliermal fermion masses are generally 
not important because for a <^ 1, T ~ M they are much smaller than M and the involved 
momenta. 

Nevertheless there are situations in which fermions with thermal masses can be relevant. 
The most obvious is F ^ g/M, h. Then the infiaton effectively only couples to fermions. 
However, due to the Pauli blocking the maximal temperature that can be achieved by 
perturbative decay into fermions is usually smaller than the infiaton mass and Tc (unless 
one choses a very large coupling Y, which is difficult to consort with a fiat effective infiaton 
potential and the slow roll conditions), thus the effect of thermal fermion masses is again 
negligible. 

Another possibility is that the effective masses for bosons grow faster with temperature 
than those of the fermions. Then the phase space for the decay into bosons closes while 
that for fermions is still open. If all couplings are weak the temperature at which this 
happens is larger than M, the previous arguments apply and the decay rate into fermions 
is small. An interesting situation can occur if the thermal boson masses Mi and M2 are 
different. Then the suppressed, but nonzero decay rate into fermions can play the same 
role as off-shell effects in scenario III): It can allow the system to pass the temperature 
region Tc < T < Tc. At temperatures T > Tc Landau damping via processes Xi4' Xj is 
efficient, see figure [181 In this case thermal fermion masses are decisive because they decide 
if r in the forbidden region is large enough to heat the plasma up to Tc. If the thermal 
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Figure 17: Fq from the trilinear interaction for ai = 0.01, a2 = 0.39 and a = 1 without 
resummation (green dotted), with thermal masses (blue solid) and with thermal masses 
and widths (red dashed). 
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Figure 18: Evolution of T (solid line) and the relative contributions of (dotted) and 
(dashed) to the energy density for parameter choice IV), g/M = \/2 ■ 10~^, ai = 0.01, 
a2 = 0.39, a = 0.1. The blue lines show the result when only the trilinear interaction 
is taken into account, the red line shows the behaviour when quartic interactions with 
h/2 = g/M are added. 
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masses of and \l/2 are different there is a temperature regime in which the plasma can 
also be heated by processes \E'j. 

A situation in which fermions with thermal masses can play an important role in per- 
turbative reheating is a ~ 1. In this case they are relevant at temperatures Tc ~ M. 
Then it is indeed possible that the decay into fermions gives a comparable contribution as 
for bosons. However, in this situation one would generally expect that the widths of the 

are not negligible. Then off-shell effects are relevant, e.g. 2 — )■ 2 scatterings involving 
intermediate off-shell \E'j. 

The efficiency of processes involving fermions can be increased by going beyond the 
quasiparticle regime. It is easy to see that cuts through the self energy insertion in the 
diagram shown in figure [10] can leave pieces with only bosons in the final state. Once 
transitions involving many quanta contribute significantly, the rate for processes involving 
fermions can be amplified by large occupation numbers in the bosonic sector. 

7.3 Discussion of Thermal Masses during Reheating 

Thermal masses of the decay products can modify the rate at which the inflaton releases 
energy into the primordial plasma during reheating. This can impose an upper bound on 
the temperature if several conditions are fulfilled. First, thermal masses have to dominate 
the effective masses. This is only the case if the amplitude of the (0) oscillations is not too 
large (e.g. a^T^ > for the scalar trilinear interaction). The applicability of our analysis 
is restricted to perturbative inflaton decay. 

Second, the dissipation process has to be driven by decays (or scatterings) of sin- 
gle (quasi) particles. At high temperature, contributions from processes involving many 
quanta become increasingly important. In our setup the contributions from these are 
parameterised by the widths of the resonances in the plasma. An upper bound on the tem- 
perature can only be imposed if all relevant fields are in the quasiparticle regime, i.e. the 
propagating states have a narrow width. Only then the dissipation can be pictured in terms 
of approximately energy conserving decays and scatterings of individual quasiparticles. 

Finally, assuming the above conditions are fulfilled, a blocking of the phase space for 
the decay only stops reheating if no other channels of dissipation are available. When (f) 
mainly couples via 3-vertices, dissipation at low temperatures is dominated by 0-decay. 
If at the temperature when the phase space for the last decay channel closes, some 
quasiparticle in the plasma is heavy enough to decay into a final state including 0, e.g. 
X2 't'Xii can dissipate energy into the plasma via the inverse process (pxi X2- 
The competition between these processes leads to an overall energy transfer from into 
the plasma analogue to Landau damping as long as the occupation numbers are larger 
than their equilibrium values. If there is no sufficiently heavy quasiparticle in the plasma, 
reheating indeed becomes inefficient for T > T^. 2 — )■ 2 scatterings are only possible 
at second order in the coupling (while the leading contribution to the decay comes at 
first order). Furthermore, in the regime where the decay is forbidden (Mi + M2 > M, 
M > |Mi — M2I), scatterings involve a virtual intermediate quasiparticle and are suppressed 
by the corresponding energy scale. However, when couples to other fields via 4-vertices, 
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2 — )■ 2 scatterings between real quasiparticles are always possible at leading order. They 
require real bath quanta in the initial state, thus their efficiency increases with temperature 
(they are negligible for T ^ M). The same applies to effective vertices with more legs. 
Thus, thermal masses only imply an upper bound on the reheating temperature if they 
dominate the effective masses, the propagating states in the plasma are narrow resonances 
and scatterings and Landau damping are negligible. 

These conclusions can be obtained from the results of the previous sections. We also 
demonstrated them numerically and computed the time evolution of the temperature for 
several parameter choices. We would like to emphasise again that these parameter choices 
are not meant to represent a particular (realistic) model of inflation, but were chosen for 
illustrative purpose. In models of single field inflation there exist strong bounds on the 
couplings of (p to other fields (and weaker bounds on the interactions of the decay products) 
- too strong couplings would spoil the fiat effective potential necessary for slow roll inflation 
by radiative corrections. The field (p may also be viewed as the relevant direction in field 
space during reheating in multi-field models. In this situation the bounds are weaker. 
However, our goal here is not to impose bounds on the parameters of a particular class of 
models, but to perform a self-consistent study of a potentially relevant effect. We leave the 
application to realistic scenarios to later work. 

Finally, we would like to emphasise that in reality the details of the reheating process 
are far more complicated than in our model. Even during perturbative decay there are 
various different time scales related to the equilibration of the decay products. As long as 
all of them are much shorter than l/F we do not expect these details to affect the rate at 
which (f) dissipates its energy considerably. We discuss these aspects at the end of section 
1831 

8 Discussion and Conclusions 

We studied the relaxation of scalar and fermionic fields in a large thermal bath from first 
principles of nonequilibrium quantum field theory. This situation is realised in various 
interesting physical systems, including thermal production of particles from a plasma, 
propagation in dense matter or cosmological freezeout processes. Depending on the initial 
state of the out-of-equilibrium fields, relaxation can mean dissipation of energy into the 
bath or thermal production of quanta. The formalism we used allows a quantum description 
of both processes. 

We performed all computations in terms of correlation functions for quantum fields, 
avoiding the ambiguities in Boltzmann type kinetic equations related to quantum coher- 
ence, the definition of asymptotic states or particle numbers in a dense plasma and the 
assumption of molecular chaos. Basis of our analysis are exact expressions for nonequilib- 
rium correlation functions for arbitrarily large deviations from thermal equilibrium that are 
valid for all times. They were obtained without semiclassical approximation or a gradient 
expansion. 

The main goals of our study were to explore which effect the modified dispersion rela- 
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tions in the plasma have on the relaxation rate and under which circumstances a description 
in terms of quasiparticles is suitable for a relativistic quantum system. 

8.1 Summary of Results 

At high temperature and density, relaxation in a plasma is not only driven by decays and 
scatterings of individual particles. Processes involving many quanta, that naively are of 
higher order in the coupling, have to be taken account. Their contribution is negligible in 
a dilute plasma, but grows with the occupation numbers. Therefore consistent results can 
only be obtained using resummed perturbation theory. In our setup the contribution from 
processes involving many quanta can be parametrised by the widths of the resonances in 
the plasma. 

If the widths of all resonances are small, these can be interpreted as quasiparticles with 
dispersion relations given by the poles of the resummed propagators. Then the exchange 
of energy between the modes of different fields can effectively be described in terms of 
decays and scatterings amongst individual quasiparticles. This is an essential condition 
for the validity of effective kinetic equations. However, even in the quasiparticle regime 
it is generally not valid to obtain such equations by simply replacing the vacuum masses 
by thermal masses because the dispersion relations can be more complicated. In addition, 
collective excitations may appear as propagating states. 

Though effective kinetic equations for quasiparticles describe the energy transfer be- 
tween different degrees of freedom with good accuracy, they do not capture all properties of 
the system. We explicitly showed that the energy density of the out-of-equilibrium field is 
not that of a quasiparticle gas and does not follow a Boltzmann equation for quasiparticles. 

The modified dispersion relations in a plasma lead to a number of interesting effects. 
If the out-of-equilibrium field couples via 3-vertices, the leading order contribution to 
the energy exchange with the bath can effectively be switched off in some temperature 
regime. This happens when, due to the temperature dependent dispersion relations, none 
of the quasiparticles connected by the vertex can kinematically decay into the other two. 
This can considerably delay the relaxation and have interesting cosmological implications, 
for instance during reheating or for the thermal production of particles. Interactions by 
4-vertices always allow for relaxation at leading order as these mediate 2 — 2 scatterings. 
However, since they require real quanta from the plasma in the initial state they give only 
small contributions at low temperature. 

The widths of the quasiparticles parametrise the effect of processes that involve more 
quanta and vertices. These appear to be energy violating in the quasiparticle picture, 
but indeed only reflect the fact that relaxation can also happen via off-shell processes 
and multiple scatterings with many quanta from the bath. At high temperature there 
are various effects including induced transitions for bosons, the large density of scattering 
partners or the Landau-Pomeranchuk-Migdal effect which can make these contributions 
comparable to or even larger than the rate at leading order. Mathematically this manifests 
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in the poor convergence of the perturbative series and the need for resummation. In the 
quasiparticle picture, this appears as an increase of the width or melting of the resonance 
peaks at high temperature. Physically it means that in a dense plasma the picture of 
individual particles travelling freely between isolated interactions does not hold. 

When the widths of the resonances are not small, the quasiparticle interpretation and 
description in terms of kinetic equations break down. For instance, even in a temperature 
region where all leading order processes involving real quasiparticles are kinematically 
forbidden, the relaxation rate may be as big as in the region where decays and inverse 
decays are allowed. Thus, naive estimates based on kinematics involving thermal masses 
fail to even qualitatively reproduce the behaviour. In the quasiparticle picture this can 
be interpreted in terms of off-shell processes, mediated by virtual quasiparticles. At low 
temperatures, these effects can be included into the Boltzmann approach by adding matrix 
elements for higher order processes. At large temperatures, this would require the inclu- 
sion of infinitely many diagrams with large numbers of particles in the initial and final state. 

We computed analytic expressions for the relaxation rates of scalars and fermions via 
a trilinear and a Yukawa interaction, respectively, for arbitrary energies and momenta in 
the quasiparticle regime. We furthermore derived a formula for the rate of relaxation via 
the trilinear scalar interaction, expressed as a one-dimensional integral, which includes the 
effects of higher order processes. Comparing to the quasiparticle result, we found that the 
corrections can be large even for a moderate width Fj/Mj ~ 10~^. 

We applied our results to a particular cosmological problem, namely the question 
whether large thermal masses of the decay products can close the phase space for pertur- 
bative inflaton decay and therefore impose an upper bound on the reheating temperature. 

Finally, we would like to emphasise that in reality the details of the relaxation process 
can be far more complicated than in our model. We assumed instantaneous equilibration 
of the bath degrees of freedom. There may be different time scales related to this process. 
The word "quasiparticles" has been used for propagating states with definite momentum, 
the modes of a homogeneous field, with properties defined by averaging over a statistical 
ensemble. They are not localised particles (wave packets). When the relaxation is pic- 
tured as a statistical process, driven by the competition of decays and inverse decays (and 
scatterings) of quasiparticles, each individual of these processes is localised in time and 
space. Even though the relaxation time of is much larger than that of the other fields, 
the details of the thermalisation process and local deviations from equilibrium do affect the 
decay of individual quasiparticles. However, we are not interested in the fate of individual 
quasiparticles here, but in the overall relaxation of the system. Therefore we believe that 
the thermal ensemble average is justified and gives correct results as long as the relaxation 
time in the plasma is much shorter than 1/F 
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8.2 Conclusions and Outlook 



Boltzmann equations and their matrix valued generalisations can accurately describe many 
nonequilibrium phenomena in a medium. They are known to suffer from uncertainties when 
the coherence lengths are large, but are usually assumed to be accurate in absence of such 
effects. On the other hand, it has been known for long that resummation is necessary in 
gauge theories at high temperature because processes involving many quanta that naively 
are of higher order in the coupling in fact contribute at leading order. It is important to 
understand in which cases and how these many-particle amplitudes can be incorporated 
into the framework of single particle distribution functions. 

The inconsistencies are related to the assumption that the system should be described 
in terms of individual particles as asymptotic states. We showed that they can be removed 
in an intuitive way when swapping the phase space distribution functions as dynamical 
variables for correlation functions of quantum fields. 

The Schwinger-Keldysh formalism used in this work allows for a full quantum treatment 
of nonequilibrium phenomena in the framework of quantum field theory. A consistent 
perturbation theory can be formulated in the framework of the nPl effective action. It 
has previously been applied to problems in which all degrees of freedom are far from 
equilibrium, including preheating after inflation and equilibration after heavy ion collisions. 
However, the equations of motion, coupled second order integro-differential equations, are 
considerably more complicated than the semiclassical Boltzmann equations. In most cases 
they can only be solved numerically in toy models. The use of exact equations comes at 
a price, the loss of the simplicity, transparency and intuitive physical interpretation that 
make the Boltzmann equations appealing. 

In this work we studied systems in which only few degrees of freedom are out of equilib- 
rium and weakly coupled to a large thermal bath. In this setup most of the computations 
can be performed analytically and the quantities in the field theoretical approach have an 
intuitive physical interpretation. The analytic solutions allowed us to study in detail the 
effect of quantum and higher order corrections, parametrise them in quasiparticle prop- 
erties and understand the limits of this approximation. The semiclassical description in 
terms of effective kinetic equations can be recovered in the limit of a dilute gas. 

The main simplification that allowed a widely analytical treatment lies in the negligible 
backreaction. We have neglected backreaction in two ways: by ignoring the effect of the 
out-of-equilibrium fields on the bath temperature and by dropping contributions to the self 
energies that have nonequilibrium propagators in the loops. This allows to use resummation 
techniques known from the theory in equilibrium. There are interesting systems in which 
the second assumption is not justified while the first one still holds. In leptogenesis, for 
instance, the baryon asymmetry of the universe originates from CP-violating interactions 
of heavy right handed neutrinos that are out of equilibrium with a thermal bath of SM 
fields. Dropping the small Cabibbo-Kobayashi-Maskawa-phase in the SM, their Yukawa 
couplings are the only source of CP-violation. When computing the CP-asymmetry, the 
suppression of contributions from diagrams that have the out-of-equilibrium fields in the 
loop by the large number of degrees of freedom in the bath does not apply as none of the 
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other interactions is CP- violating. Indeed, their inclusion is crucial as the generation of 
an asymmetry comes from a quantum interference between the two leading orders in the 
Yukawa couplings. 

Therefore it is important to develop techniques that allow a consistent evaluation of 
higher order diagrams in which the nonequilibrium fields appear in the loop while making 
maximal use of the simplifications due to the weak coupling to a thermal bath and negli- 
gible backreaction on the temperature. The nonequilibrium propagators used in this work 
can be used for a perturbative treatment and the formulation of Feynman rules is straight- 
forward. However, practical evaluation of higher order diagrams is technically difficult due 
to the dependence on centre of mass time, and it remains to be seen how to perform the 
resummation of all stronger (gauge) interactions in the bath. 
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A Appendix 



A.l Nonequilibrium Quantum Field Theory 

Here we summarise basic concepts of nonequilibrium quantum field theory and notations 
used in this work. For detailed reviews on the topic, we refer the reader to [^H3]. Most of 
the notation we use here is adopted from j6], which also contains a slightly more detailed 
introduction into the topic than presented here. 

Out of equilibrium there are two independent two point functions for each field. For a 
real scalar field these can be chosen to be the connected Wightman functions 

A>(xi,a:2) = (0(xi)</.(x2))c, (A.l) 
A<(Xi,X2) = (0(x2)0(xi)),. (A.2) 

where the (. . .) is to be understood in the sense of (12. ip and includes averaging over 
thermal as well as quantum fluctuations. The subscript c indicates that are connected 
correlation functions. For a Dirac fermion the corresponding definitions are 

S>^iXi,X2) = {^aiXi)^^{x2))c (A.3) 
5<^(Xi,X2) = -{^(3ix2)^aiXi))c (A.4) 

Here \1/ is a Dirac spinor and a and /3 are spinor indices which we will always suppress in 
the following. It is convenient to define the linear combinations 

A-(X1,X2) = 2(A>(X1,X2) - A<(X1,X2)) (A.5) 

A+(xi,X2) = i(A>(xi,X2) + A<(xi,X2)) (A.6) 
S~{xi,X2) = i{S^{xi,X2) - S'^{xi,X2)) (A.7) 

S+{Xi,X2) = 1{S>{XUX2) + S<{XUX2)). (A.8) 



The functions (lA.SP - OA.SP have intuitive physical interpretations. The spectral functions 



A~ and S~ are the Fourier transforms of the spectral densities and, roughly speaking, 
characterise the spectrum of excitations in the plasma. The statistical propagators A+ and 
can be viewed as a measure for the occupation numbers. They fulfil the symmetry 
relations 

A~(a;2,a;i) = -A~(xi,X2), (A.9) 
A+(x2,Xi) = A+(a:i,X2), (A.IO) 

S~{x2,xi) = -7o (5'~(xi,X2))So (A. 11) 

S+{x2,xi) = 7o(5+(xi,X2))So, (A.12) 
which can be seen from the definitions, and the boundary conditions 

A-{xi,X2)\t,=t, = 0, (A.13) 
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dt,A~{xi,X2)\t,=t2 = -dt,A"{xuX2)\t,=t,=S^^\^i-^2), (A.14) 
9i,9i,A-(xi,X2)|t,=t, = 0, (A.15) 
S~{xi,X2)\t^=t2 = i7o5(xi-X2), (A. 16) 

which follow from the canonical commutation and anticommutation relations. The bound- 
ary conditions for the statistical propagators are determined by the physical initial condi- 
tions in which the system is prepared. 

In thermal equilibriurrn or vacuunQ the correlation functions only depend on the 
difference of coordinates Xi — X2- Furthermore they are related by the Kubo-Martin- 
Schwinger (KMS) relations, which are most conveniently written in terms of their Fourier 
transforms in Xi — X2- In absence of chemical potentials they read 

A<e.(^) = e-^-A>,^(a;) , S<,^{u) = -e-^-S^^^iu) . (A.17) 

Here (w, q) is a momentum space four vector and /3 = 1/T the inverse temperaturj^. 
The KMS relations are a consequence of the fact that the equilibrium density matrix 
g'^'^ = exp(— /3H)/Tr exp(— /3H), where H is the Hamiltonian, acts as a time evolution 
operator in imaginary time and imposes boundary conditions on the in the complex 
time plane. They allow to express all two-point functions in terms of the spectral densities, 



A<e,(w) = /B(^)Pq(c^) , Aq,e,(^) = (1 + /b (w) )Pq(c^) , 
S^eqi^) = (1 - /F(u;))Pq(u;) , S<,^iu) = -Mu)p^{u). 



(A.18) 



(A.19) 



Here /b = {e^'^ — 1)"^ and fp = {^^^ + 1)"^- The Bose-Einstein and Fermi-Dirac distri- 
butions and fp naturally arise as a consequences of the boundary conditions for the 
correlation functions. The concrete shape of the spectral densities Pq(co') for scalars and 
Pq(a;) for fermions will be specified below. 

Out of equilibrium ( I2.2p -( l5^ all independently depend on both arguments Xi and X2- 
They also cannot be related to each other by a relation analogue to ( ]A.17P -( fA.19l) . The 



spectral function and statistical propagator both have to be found as solutions to a set of 
coupled integro-differential equations known as Kadanoff-Baym equations. 

They can be derived in the Schwinger-Keldysh formalism 1751. For simplicity we restrict 
ourselves to cases where the latter are Gaussian at initial timcll, i.e. all n-point functions 
at initial time can be expressed by the two point function. Higher order correlations 



^"^See |42) for a detailed introduction to quantum field theory in equilibrium. 

^^The correlation functions in vacuum can be obtained as the zero temperature limit of those in equi- 
librium and we in the following understand vacuum as a special case of thermal equilibrium. 

^^Here and in the following we work in the restframe of the bath, where T has a physical interpretation 
as temperature. Due to the choice of frame the expressions are not Lorentz- invariant. This choice is for 
convenience, the theory can of course be formulated in a covariant manner |25j . 

^^For non-Gaussian initial conditions the formalism described here remains applicable, but some addi- 
tional terms appear in the Kadanoff-Baym equations at initial time [50| . 
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Figure 19: The Keldysh contour runs from some initial time x° = ti + ie parallel to the real 
axis = t + ie) up to some final time tf + ie and returns to ti — ie. To compute physical 
correlation functions for arbitrary times t > ti, one takes the limits t/ — )■ oo and e — )■ 0. 



unavoidably build up at later time when system and bath are brought into touch, which is 
perceived as an increase of entropy by observers looking at Gaussian correlators only |49] . 
Starting point is the generating functional for correlation functions of fields with time 
arguments on a contour in the complex time plane known as Keldysh contour, (cf. figure 
[T9|) . The choice of this contour (rather than the real axis) is a consequence of the fact that 
nonequilibrium processes are initial value problems. The system is prepared at an initial 
time ti and its state at later times is unknown. The Keldysh-contour, which starts and 
ends at the same time, allows to define a generating functional without knowledge of the 
state at asymptotic times t = ±00. This is in contrast to the 'in-out' formalism used to 
compute the S-matrij@- For the time ordered two-point function Ac{xi, X2) for scalars on 
the contour C it yields the Dyson- Schwinger equation 

(□1 + m^)Ac{xi, X2) + J d'^x'Ucixi, x')Ac{x', X2) = -iSc{xi - X2) . (A.20) 

Here Di = d'^dfj,, where the derivatives are with respect to the components of Xi. Ilc{xi, X2) 
is the self energy and Sc{xi — X2) a 5-function on the contour. Green's function and self 
energy can be decomposed as 

Ac(xi, X2) = edxl, x°)A>(xi, X2) + 9c{xl x?)A<(xi, X2) (A.21) 

nc(xi, X2) = Ocixl, x°)n>(xi, X2) + edxl, x?)n<(xi, xa) (A.22) 



In (IA.20p the time coordinates of Ac and He can lie on the upper or the lower branch of 



the contour. When both time arguments lie on the upper branch, Ac can be identified 
with the usual Feynman propagator Aj? in the limit e — )■ 0, cf. figure [121 For both time 
arguments on the lower branch this limit yields the anti-time ordered propagator A^ due 
to the inverse direction of the ordering along the contour. In the cases that Xi lies on 
the upper and on the lower branch and vice versa, Ac can be identified with A^ and 
A^, respectively. In a perturbative expansion of flA.20p in terms of Feynman diagrams. 



the time arguments of internal vertices can lie on either branch. Hence, the number of 
contributing graphs doubles with each internal vertex since all possible combinations are 
summed ovecj. Two upper vertices are connected by Ap, two lower vertices by Ap and 



^^Due to this fact this formaHsm is sometimes called 'in-in' formalism. 
■*^This fact is sometimes referred to as 'doubling of degrees of freedom'. 
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vertices of different type by and A^. Eacli lower vertex leads to an additional factor 
— 1. This is reminiscent of the real time formalism used in equilibrium. All correlation 
functions can be obtained as linear combinations of A~ and A"*" via 
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(A.26) 



Using the above relations, a straightforward calculation allows to derive the equations of 
motion for A~ and A'^ from (lA.20p . known as Kadanoff-Baym equations, 

(□1 + m2)A"(xi,a;2) = - j d^ii! j^' dt'ir{xi,x')A~{x\x2) , (A.27) 

(□1 +m2)A+(xi,X2) = - J d^x.' J^^ c/t'n~(xi,x')A+(x',X2) 

+ [ rfV dt'U+{xux')A-{x',X2) . (A.28) 
J Jti 

Here the self energies 11 are defined analogue to 

n+(xi,X2) = -^(n>(xi,a;2) + n<(xi,a;2)) (A.29) 
n-(xi,a;2) = n>(xi,X2) -n<(xi,X2), (A.30) 
and can be related to the retarded and advanced self energies by 

n^(xi, X2) = e{h - t2)n-(xi, X2) , n^(xi, xs) = -e{t2 - ti)n-(xi, xs). (A.31) 

All corresponding equations for the fermionic Green's functions are obtained by the re- 
placements (Di + m?) -{i^i - m), Ac(xi,X2) S'c(xi,X2), nc(xi,X2) Ec(xi,X2) 
and so on. 

The Kadanoff-Baym equations (IA.27P and (1A.28P are the exact equations of motion for 
the correlation functions. They do not rely on a gradient expansion or semiclassical ap- 
proximation and make not restrictions on the size of the initial deviation from equilibrium. 
Roughly speaking. A"*" and are the field theoretical generalisations of the phase space 
distribution functions appearing in the semiclassical Boltzmann equations while the self 
energies take the role of the collision terms. This interpretation is intuitive as the total 
cross section is related to the discontinuity of the self energy by the optical theorem. 



Depending on whether one defines the self energies as the sum of all IPI-diagrams, c.f. [M], or i times 
that sum, c.f. [6], these definitions are either completely analogue to (jA.sp . (jA.6P or vary by a factor i. 
Here chose the 'asymmetric' convention used in [6]. 
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We are mainly interested in cosmological applications and restrict the analysis to spa- 
cially homogeneous and isotropic systems, in which the correlation functions can depend 
on the difference of spacial coordinates Xi — X2 only. Furthermore, we are interested in the 
case that or \E' are weakly coupled to a large thermal bath. Formally this corresponds 
to a density matrix of the form ^^ys ® ^^Ith' where gH^^ = exp(-/3Hbath)/Tr exp(-/3Hbath) 
and Qsys is the initial density matrix of the system. The interactions of the fields that 
make up the bath amongst each other are much stronger than the couplings to or \& 
and bring them to equilibrium fast on the time scale associated to the dynamics of the 
out-of-equilibrium degrees of freedom. Then all correlation functions of the bath fields can 
always be characterised by a single temperature T (and possibly a handful of chemical 
potentials). We are interested in the case that the out-of-equilibrium fields mainly interact 
with (in particular decay into) bath fields only. This can explicitly be realised by using 
coupling terms that are linear in (p and \E' and can be written in the form a(f)0(p[X] and 
Q!^'0^[A'] + h.c Here 0^\X\ and are composite operators that contain bath fields 

only, to which we collectively refer as and a, a are coupling constants. The effects 
of backreaction on the self energies are suppressed by additional powers of the coupling 
constants as well as the large number of degrees of freedom in the bath and can safely be 
neglected. Thus, the self energies are determined by bath properties. In particular, they 
depend on the relative coordinate x\ — X2 only, inherit the KMS relations 



n>(a;) = e'3-n<(a;) S>(a;) = -ef^-J:<{u) 

and can be computed perturbatively by methods of thermal field theory. It can be proven 
that in this case also the spectral functions are time translation invariant [6]. 
The definition (lA.3ip gives rise to the well-known spectral representation 

n^(c.) = ^ r ^^^n^ (A.33) 



which implies 

nV) = in-M + p/|i^Il;M (A.34) 

and an analogue relation for fermions, where V denotes the principal value. These relate 
the discontinuities of the self energies to 11" and S~. For a real scalar with real couplings 
fEOil) impliesEl] 

-Imn^(a;) = ^n-(a;) . (A.35) 

The quantities 11^ individually can be related to the gain and loss terms in fl2.25p . which 
fulfil the detailed balance ratio due to (1A.32P when the bath is in equilibrium. 



"^A table of the analytic properties of the self energies can be found in the appendix of [B]. 
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In the mixed time-momentum representation the Kadanoff-Baym equations can be 
written in the convenient form 

{dl+ujl)^-{h-h) = - dt'n-(ti-t')A-(t'-t2) , (A.36) 

Jti 

+ I ' c/t'n+(ti - t')A-{t' - t2) . (A.37) 
Ju 



{ilodt,-<^-xa)S^{h-t2) = r dt'T,-^{h-t')S^{t' -h) , (A.38) 

Jt2 

Jti 

- dt'mt,-t')S^{t'-t2) . (A.39) 

where cjq = (q^ + m^)^/^, ^ = 'jiQi and ti is the initial time at which the system is prepared, 
which can be set to zero without loss of generality. 

All formulae presented here are valid in Minkowski spacetime. They are valid in good 
approximation if the characteristic time scale of the nonequilibrium dynamics is short in 
comparison to that of Hubble expansion. A generalisation to curved spacetime has e.g. 
been discussed in 



A. 2 The Fermion Self Energy 

For vanishing m (I5.15P reads 

d^p 1 



(27r)2 2(]i2r22 

1 - MQ^) + /b((^2)) [{^iio -in)K^ -^i- ^2) 

+ (/^(fii) + /b(1^2)) ((1^170 -|)7)5(u; - 1^1 + ^2) 

+ (l]i7o + pj)5{lo + n^-n2))] (A.40) 



We take a; > 0, which is sufficient due to the symmetry of fl5.12p . The self energy can be 
decomposed as Sq(a;) = 2i{aq{u)(/+ bq{uj)i/), where u is the four velocity of the bath and 
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aq(a;), &q(w) are scalar functions that depend on u and q. In the restframe of the bath 
u = (1, 0). We define the projections 

M^) = ^tr(gS-(a;)) (A.41) 

B^iuj) = ^tr(^q(c^)) (A.42) 
which can be related to the coefficients aq(a;) and &q(u;) via 

[quy — q^u^ q"' 

2a,(.) = + t'^")"" ^ "■^-<"' - "'^-<") + (A.44) 
[qu)'^ — q^w^ q^ 

From this one can find 

+ (wfil + qp)6{uj + 1^1 + i]2)) 
+ + /2) ((c^f^i - qp)5(u; ^2) 

+ (c^(]i + qp)5(a; + fii-fi2))) (A.45) 

with /i = and /2 = fsiS^i) during this computation. This expression can be 

rewritten as 



A = -Y^ 



-f\ + h){S, + 52) + (/i + f2){Ss + 54)) 
-qp((l -fi + /2)(5i - ^2) + (/i + /2)(54 - 53))) (A.46) 

where 

5^ = - f2i - ^^2), ^2 = 5(w + i]i + ^2), fA 47'! 

53 = 5(^ + 1^1-1^2), 54 = 5(^-^1 + 1^2) ■ 

We drop 52 as it cannot contribute for a; > and change to spherical coordinates ip, |p|. 
The if integration is trivial. As the fermion in the loop is massless |p| = Cli. Introducing 
X = |p||q| cos{'&) = pq one can write 

A = — - / dCli / dx 

16vr|q| Jo J-n,\ci\ 

( 6{x-Xol){{uJn^-x){l-f^ + h)) 
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= liq'-co' + M') 
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— Qiu 


Xo4 







(A.48) 



(A.49) 

(A.50) 
(A.51) 

This allows to perform the x integration, 

A = ( ^ dn, {tufi j^d^i {cofs - ^3) + ^ {ufi + ^4) ) ■ (A.52) 

Here the subscript at the J. indicates which Si determines limits for the Qi integration. 
The fi and ji are given by 

(A.53) 
(A.54) 
(A.55) 

(A.56) 
(A.57) 
(A.58) 

It is easy to see from Eqs. ([OT]) . (1X42|) and (1X48|) that 

B = (^J^dn,f, + dn,f, + ^ ^^^1/4) (A.59) 

/i and £1 are contributions from decay and inverse decays of \E' and can lead to a zero 
temperature part in case (a), /a, /4, £3 and ^4 are related to Landau damping. Note that 

/i = -/4 (A.60) 
£1 = (A.61) 

due to the property fsi—z) = — 1 — fsiz) despite the fact that the terms originate from 
different processe^j. The fi are symmetric in u, the £i antisymmetric. B is symmetric in 
u} while A is antisymmetric. As a consequence, a is antisymmetric while b is symmetric. 
The stem functions of all /j, £i are known analytically (see flA.75P - (lA.80|) ) and the only 
remaining task is the determination of the integration limits for each 5j-function. 



/l 






/s 






/4 












£3 






Si 


= {^{u^ - q2 - M^) - {MQi) + fB{ni - 


-UJ)) 



^^Note that despite the equahties (jA.60p and (|A.6ip the Landau damping terms /j, ^4 never lead to a 
contribution to S at zero temperature while the decay and inverse decay parts /i and ^1 can contribute 
as expected. The reason lies in the different integration limits, see (|A:62|1 and (|A7T|l 
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6{u — Qi — Q2) : In order for xqi to be a zero, the condition 

w - fii > (A.62) 

has to be fulfilled. In any case, 

1^1 > m = 0. (A.63) 
In order for the x-integral to be non-zero it is required that 

|xoi| < ^i\q\. (A. 64) 

The solutions to |xoi| = fii|q| are 

d;± = -^— ^(a;±|q|) (A.65) 
z q 

with = cj^ — q^. One has to distinguish three different regimes: For < w < |q| and 
f2i > only Cj+ is a solution, and it imposes a lower bound on Qi in order for the inequality 
([Oil) to be fulfilled, leading to > a)+. On the other hand the condition flA.62p has to 
be fulfilled, and since for u < |q| always uj^ > u, there is no contribution to the integral 
from this region. For |q| < u; < (q^ + M^) 2 both u}± < and none of them makes (1A.69P 



an equality. For u > (q^ + M^)2 both u± are always smaller than u and (1A.69P leads to 
u!^ < uj < u+. Therefore 

[ d^i = e{q^ - M^) f ^ d^i (A.66) 

6{uj + Qi — Q2) '■ Here the three conditions 

(A.67) 

m = (A.68) 

fii|q| (A.69) 





> 




> 




< 



have to be fulfilled. In this case (1A.69P is made an equality for Qi = —u}±. The same 



regimes have to be distinguished. For a; < |q| only —U- makes (1A.69P an equality while 
— a)_i_ is negative and not a solution. — a)_ is positive as required by (1A.68P and forms a 
lower bound. For |q| < < (q^ + M^)2 both —Cj± are positive and solutions. Beyond 
its first order pole at u = |q| the solution — a)+ is now the larger one and forms an upper 
limit, leading to — cl;_ < Qi < — a)_|_. For u > (q^ + M^)2 both — w± are negative and not 
solutions of ( ]A.69P as an equality. Then there is no contribution to the integral from that 



region. Therefore 



/ rffii = 9{-q^) / c/fii + 9{q^)9{AP - g2) / (A.70) 

J 3 J —uj— J —Hi— 
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6{u — Qi + Q2) '■ The situation here is exactly the same as for 61, in particular Xoa = Xqi, 
except that the condition u — Qi > has to be replaced by 



w - fii < 0, 



(A.71) 



enforcing Cli > u. Again for u < |q| only a)+ fulfils (lA.64p . imposing a lower bound on Qi 
and for u > (q^ + M^)2 both oj± are solutions. a)+ is the upper and w_ the lower bound 
here. For < < none of u}±. is a valid solution. This time the condition (lA.7ip 
selects out the region < 0, hence the integral is 



U1+ 
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(A.72) 



(A.73) 



(A.74) 



(A.75) 



(A.76) 
(A.77) 



J' 



Li2(e^(^^--)) -Li2(-e^''^)) (A.78) 



g, = ^^^^(in (1 + e^^^) - In (e^(^^+-) - l) 



2(3 

^4 = -gi. (A.80) 

Li2 is the dilogarithm function. The jTj as displayed here are not real in general due to the 
choice of different branches of the (di) logarithms, but the imaginary terms always cancel 
since the choice of branch is always the same at both integration limits. This analytic 
result for S~ is in agreement with numerical plots shown in [66]. The ^-functions are, as 
in (13. 6p . a consequence of energy conservation. 
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